d) Instituto Nacional de Meteorologia e Geofi- 
sica: recolhe informação meteorológica e 
possui uma rede de radiação solar que abran- 
ge todo o país. 

e) Direcção-Geral dos Combustíveis, possui um 
grupo de trabalho (G. T. E. S.) que vem 
realizando um louvável esforço para concre- 
tizar colaborações de organismos vários no 
uso da energia solar. Tem ensaios realizados 
sobre colectores planos nacionais e estran- 
geiros, em condições reais de funcionamento. 

f) Direcção-Geral das Construções Escolares: 
tem algumas instalações em funcionamento 
e experiência de utilização em escolas, 

9) Escola Superior de Belas Artes do Porto: 
iniciou com êxito a concepção e integração 
de realizações tecnologicamente simples na 
viabilização de novas formas de «habitat». 

h) Faculdade de Economia da Universidade do 
Porto: criou grupos de trabalho que se de- 
bruçam sobre o impacto na região dos dife- 
rentes modos de conceber o desenvolvimento 
tendo em conta as potencialidades nacionais. 
As alternativas energéticas são um dos te- 
mas. 

à) Instituto Universitário de Evora: possui um 
Departamento de planeamento biofísico e ini- 
ciou estudos base relativos às potencialidades 
locais de utilização de energia solar. 


A listagem e a caracterização não são de modo 
algum exaustivas e são apenas as que naturalmente 
me ocorrem devido a conhecimento pessoal da sua 
existência, 


Apesar do seu carácter incompleto, dão já uma 
ideia das entidades oficiais que manifestaram inte- 
resse concreto avançando para estudos ou realizações 
preliminares. Apesar disso, o mercado nacional das 
instalações para aquecimento solar é ocupado na sua 
quase totalidade por equipamento importado... e na 
maioria dos casos mal adaptado e ainda pior assis- 
tido. Os importadores sabem vender, mas a sua capa- 
cidade de resposta a questões técnicas precisas quan- 
to ao equipamento é praticamente nula. O crescimento 
do mercado tende a modificar a situação. Tal como o 
aumento de taxas na importaçãão e alguma concor- 
rência que já começa a surgir tenha levado a alguma 
subcontratação de componentes à indústria nacional... 
e a descida de preços. 


Analisando um pouco mais o porquê de tal situa- 
ção num país reconhecidamente em crise ocorre per- 
guntar: é tecnologia de ponta ou avançada técnica 
de venda? A resposta mais adequada será talvez a 
de que nem uma coisa nem outra, É apenas mais 
uma manifestação entre muitas outras do trágico 
desajustamento entre a capacidade inovadora de mui- 
tos portugueses e a descrença nela posta pelos que 
dispõem dos meios de decisão que a torna operativa. 
No fosso entre ambos está a brecha por onde entra 
o produto estrangeiro... porque é estrangeiro... por- 
que se falhar o produto mas o produto for estrangei- 
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ro, ou de concepção estrangeira a desculpa é fácil... 
não se podia fazer melhor... era o que de melhor 
o estrangeiro tinha! 

No aproveitamento da energia solar para aqueci- 
mento, a incapacidade para digerir as linhas de evo- 
lução neste domínio, leva a aguardar que no estran- 
geiro se imponha... para depois se começar a pensar 
em utilizá-lo... de origem externa! 

Entretanto, uns quantos pioneiros insistiram. Uns 
quantos nacionais compraram por conta própria. Os 
jornais vão dando a notícia... e alguns (poucos) res- 
ponsáveis descobrem, subitamente e estupefactos que 
numa habitação vizinha, que no Algarve, na Praia 
das Maçãs, na Espanha onde estiveram, ou no país 
donde chegam, que aquilo no alto do telhado, era 
um colector solar e aquela esquisita arquitectura onde 
havia placas negras... se inspirava afinal num reen- 
contro com o Sol! 

Que o traço forte não choque e a ironia não fira. 
Pretendo apenas que a doutorisse acabe e se evolua 
para a competência... que é sobrevivência e pro- 
gresso. Erros assumidos, são erros ultrapassados. 


7. CAPACIDADES IMEDIATAS DE ACTUAÇÃO 


Colectores solares e termoacumuladores norma- 
lizados encontram-se já no mercado português. O ven- 
dedor estrangeiro está nalguns casos a importar ape- 
nas a superfície colectora... que naturalmente sobre- 
valoriza. Esta normalização facilita a venda, e faci- 
lita a correcção adaptativa no utilizador. Se uma uni- 
dade for insuficiente, junta-se outra... e outra. 


Encaremos esta via. É a via da produção em 
série. Qualquer metalomecânica com capacidade de 
corte e dobragem de chapa, dispõe de equipamento 
para a caixa. Para a superfície colectora é possível, 
na maioria dos casos, estudar a forma mais adequada 
ao equipamento de que se dispõe. O rendimento dessa 
superfície, seja pintura, seja tratamento químico, seja 
deposição sob vácuo, possui suficiente latitude para 
que se encontre a solução mais adequada aos meios 
fabris sem que seja necessário mais investimento 
significativo em equipamento. A eficiência térmica 
final variará, naturalmente, com a solução escolhida. 
Mas a viabilidade económica não se exprime pela 
eficiência térmica final mas sim pelo custo final da 
energia aproveitada, tendo em conta a durabilidade 
do equipamento. 


Referimos a solução metálica para a caixa. Mas 
nem a caixa nem o colector têm necessariamente que 
ser metálicos. Outros materiais poderão alcançar o 
mesmo efeito, Trata-se de encontrar a melhor solução 
entre o disponível no mercado. 

O colector normal exige vidro em chapa, ou co- 
bertura transparente com características ópticas nem 
sequer muito exigentes. Temos produção nacional de 
chapa de vidro e a generalização do aquecimento 
solar, que pode representar dezenas de milhar de 
metros quadrados por ano é um mercado importante 
para aquele sector industrial, 
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A instalação exige isolamento térmico. Temos 
isolantes nacionais e isolantes importados. Estes serão 
eventualmente, e para já, a solução mais expedita. 
Mas a incorporação em larga escala ou na totalidade 
de isolantes nacionais obtém-se em curto prazo com 
uma adequação do projecto. Os restantes acessórios 
existem no mercado e são nacionais: são canalizações, 
etc. A gama de soluções possíveis é tal, mesmo ao 
nível duma produção em série, ou semi-série, que é 
virtualmente impossível listá-las. 

Para cada caso, optimizar a solução pelos pro- 
cessos tradicionais pode ser enfadonho e moroso, Foi 
por isso que no €. T. A. M.F.U, L. se automatizou 
já o processo para os colectores. O computador existe 
para isso: para optimizar a solução e... poupar tempo 
e trabalho. 

Consideremos seguidamente o caso da adaptação 
do aquecimento a grandes imóveis ou equipamentos 
colectivos. Pode, e muitas vezes se justificará pela 
sua importância e custo, num estudo particular; tal 
não impede que muitos aspectos parcelares sejam 
normalizados. Trata-se, sobretudo, de estabelecer 
uma metodologia de projecto, de juntar um con- 
junto de informações acerca das características dos 
materiais disponíveis, capacidades de produção e pra- 
zos de entrega. Temos empresas e equipas capazes 
de fazer o projecto de águas quentes, etc. Basta que 
incorporem nos seus conhecimentos alguns aspectos 
fundamentais, que criem experiência começando, que 
agreguem arquitectos, que criem uma Tilosofia de con- 
junto, interdisciplinar... e que multipliquem o esforço 
investido automatizando muito cálculo enfadonho que 
os desanima. 

Para todos, na indústria que fabrica, na empresa 
que projecta, existem para lhes dar apoio nos casos 
singulares, para os manter a par das inovações, para 
lhes resolver as situações que exigem mestria na ciên- 
cia fundamental, as instituições universitárias. Estas, 
por seu turno, desejam que lhes transmitam a expe- 
riência que só a realização concreta das obras faculta 
e que lhes levem os problemas que o trabalho do dia 
a dia não permite aprofundar. Numa palavra, preci- 
samos de um esforço conjunto! Precisamos que não 
haja feudos, falsos prestígios, pseudocompetências. 

Capacidades imediatas existem e em larga escala, 

Tratando-se da introdução a um debate, creio 
que a sucinta exposição anterior contempla os aspec- 
tos que comecei a afirmar não deverem ser esqueci- 
das na actual situação do país. No entanto, desejaria 
sintetizar as respostas: 


Um programa de âmbito nacional de aproveita- 
mento da energia solar pode contribuir significativa- 
mente para: 


— a redução da dependência energética, 

— o equilíbrio da balança de pagamentos no tri- 
plo aspecto da redução no consumo de com- 
bustíveis, do equipamento importado e da ex- 
portação de equipamento de fabrico nacional, 

— à criação de um número elevado de empregos 
na indústria e nos serviços, 
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— à melhor utilização de equipamento improdu- 
tivo ou sub-aproveitado, 

— à melhoria da habitação social, 

— levar os benefícios da utilização de energia a 
imensas áreas do interior e estimular aí uma 
actividade produtiva na mecânica ligeira, 

— aumentar a produtividade agrícola. 


8. PROGRAMA DE ACTUAÇÃO 


Recapitulemos primeiro alguns dados actuais so- 
bre a situação portuguesa: 


a) Baixíssimo grau de utilização de energia so- 
lar apesar do enorme potencial de que dis- 
pomos. 

b) Conhecimento científico e técnico nas áreas 
críticas de utilização generalizada da energia 
solar. 

c) Equipamento fabril subaproveitado ou não 
utilizado capaz de cobrir todos os sectores 
da fabricação de equipamento. 

d) Enorme mercado potencial nos equipamentos 
colectivos, na habitação particular, na indús- 
tria e na agricultura (superior a 2 — 3 mi- 
lhões de contos ano). 

e) Crescente penetração de equipamento estran- 
geiro mal adaptado ao país. 

f) Péssimo aproveitamento das características 
térmicas naturais dos edifícios e inadequação 
destes ao clima. Pouca ou nenhuma atenção 
ao isolamento e à qualidade dos materiais 
utilizados. 

9) Ausência duma visão integrada dos proble- 
mas da energia em Portugal. 

h) Enorme dispersão de esforços, de sobreposi- 
ção de trabalhos e de equipamento. Prolife- 
ração de falsos prestígios e actuações de fa- 
chada. Lacunas de que ninguém cura mas 
são fulcrais. 

i) [Enorme desinteresse oficial pelo assunto ao 
nível dos escalões decisivos de actuação que 
se reflecte no fraco apoio aos grupos que já 
revelaram capacidade científica e técnica 
para encontrar soluções e as passarem à prá- 
tica. 


Como um programa de acção no domínio da ener- 
gia solar deve ser necessariamente encarado no âm- 
bito de um programa geral de energia; como certas 
formas de energia têm em Portugal uma situação 
previlegiada com organismos, feudos e serviços autó- 
nomos que deliberadamente, ou por simples deforma- 
ção tendem a perpetuar e aprofundar o desequilíbrio, 
a primeira medida concreta e a muito curto prazo é: 

— Criação de um Instituto Nacional de Energia 
que agrupe cientistas, técnicos, equipamentos e insta- 
lações dispersas nos Ministérios da Indústria e Tecno- 
logia (L, F. E. Nucleares, D. G. €, I. N. I. I, etc.), 
da Agricultura e Pescas, do Comércio, da Educação, 
etc. 
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Teriam participação nos orgãos de decisão as 
organizações profissionais e as empresas, etc. 

A função deste Instituto seria a de investigação 
aplicada, formação profissional a todos os níveis, di- 
vulgação, dinamização e implementação de soluções 
práticas. Não teria capacidade de imposição sobre 
as actividades económicas nem poderes coercivos. 
A sua finalidade seria de promoção tecnológica e de 
integração de esforços; não de fiscalização burocrá- 
tica. Teria algumas atribuições análogas às do L. N. 
E. €C. mas iria muito mais além nas obrigações para 
com a comunidade. A base fundamental do seu finan- 
ciamento seria uma percentagem sobre os benefícios 
económicos directos que promovesse, 

No momento actual, a criação de tal Instituto 
não aumentaria os encargos financeiros do Estado. 
Pelo contrário, a melhor gestão e dinamização dos 
recursos humanos e materiais disponíveis, permitiria 
economias e uma eficácia de actuação que a pulve- 
rização e o oportunismo não possibilitam, pelo desã- 
nimo que provocam nos mais capazes, 

Se a razão fundamental do actual estado de coi- 
sas é a ausência de uma clara política energética e a 
incapacidade para uma efectiva e incisiva actuação 
que extirpe o feudalismo e a pulverização de esforços; 
se um Instituto Nacional de Energia pode contribuir 
significativamente para ultrapassar tal situação, a 
sua existência não é indispensável para que se inicie 
desde já uma actuação concreta a qual só teria êxito 
congregando todos os esforços. Realisticamente, é pre- 
ferível ter em conta que a simples possibilidade de 
criar um Instituto Nacional é, infelizmente, susceptí- 
vel de desencadear movimentos e choques de influên- 
cias mais interessados em pelouros e poder do que em 
promover a recuperação económica e a capacidade 
tecnológica nacional. 

Feito este parêntesis, consideremos um plano pos- 
sível: 


A — ACTUAÇÃO IMEDIATA (— 6 meses) 


1 —Criar uma comissão coordenadora com re- 
presentantes das empresas de engenharia, 
do sindicato dos engenheiros técnicos e dos 
engenheiros, dos arquitectos, da indústria 
ligeira, e dos centros oficiais que se têm 
dedicado ao assunto. 

Criar uma infraestrutura de apoio nu- 
ma empresa de engenharia (de preferência 
nacionalizada) que assegure a eficácia das 
acções. 

Pela sua posição de charneira, o inte- 
resse das empresas de engenharia é fulcral. 

2 — Centrar esforços nos colectores planos. 

3 — Lançar cursos de formação e reciclagem, 
utilizando todos os meios disponíveis nas 
Escolas Superiores e Técnicas para divulgar 
o conhecimento existente sobre o assunto 
susceptível de aplicação imediata não ape- 
nas mo domínio da energia solar mas da 
energia em geral. 
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4 — Sensibilizar os grandes utilizadores (Habi- 
tação, Saúde, Ensino...) para a importância 
do balanço térmico das suas instalações 
actuais e futuras e para as potencialidades 
da energia solar. Conseguir a garantia da 
sua progressiva utilização como factor 
essencial de criação de um mercado interno 
impulsionador de toda uma nova actividade 
económica. 

5 — Iniciar a realização duma instalação protó- 
tipo de pequena escala, susceptível de gran- 
de produção standardizada. 

Iniciar o estudo e projecto duma insta- 
lação tipo, para demonstração, vulgarização 
e estudo de novas soluções. O local mais 
adequado é uma Escola de Engenharia. 

Iniciar o estudo de um projecto inte- 
grado numa instalação de grandes dimen- 
sões. 

Solicitar o financiamento através dos 
fundos para poupança de energia, de apoio 
às P.M. E,, etc. 

6 — Sistematizar e divulgar toda a informação 
já existente sobre os aspectos climatológi- 
cos do país relevantes para a utilização ra- 
cional da energia solar. 

7 — Inventariar as capacidades produtivas da 
indústria. Estabelecera padrões de qualida- 
de, Solicitar aos fabricantes características 
térmicas e mecânicas dos seus produtos. 
Conseguir da indústria da chapa de vidro a 
determinação e manutenção das caracterís- 
ticas ópticas. 

8 — Examinar, discutir e apreciar as realizações 
já existentes no país. Aferir objectivamente 
a sua adequação e integrar tal conhecimen- 
to em projectos a lançar. 


B — CURTO PRAZO (— até 2 anos) 


1 — Estabelecer claramente as potencialidades 
de cada fabricante e os padrões de quali- 
dade que é capaz e deseja manter, 

2 — Prospectar o mercado externo. Determinar 
as características mais desejáveis do equi- 
pamento. 

3 — Optimizar projectos que incorporem produ- 
tos nacionais. 

4 — Estabelecer normas de qualidade e segu- 
rança. 

9 — Estabelecer uma metodologia tipo de pro- 
jecto para pequenas e grandes instalações. 

Automatizar tudo quanto é repetitivo. 
Usar as enormíssimas capacidades disponí- 
veis no parque nacional de computadores. 
Estimular a permanente actualização das 
metodologias e métodos recorrendo aos or- 
ganismos de investigação especializados. 

6 — Criar um banco de dados sobre materiais, 
características e preços. Usá-lo como ins- 
trumento de análise de importações. 
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7 —Ensaiar, testar e desenvolver sistematica- 
mente movos protótipos, Concretizar indus- 
trialmente as soluções mais promissoras en- 
contradas nos laboratórios de investigação. 

8 — Desenvolver as aplicações na refrigeração e 
articular a acção com o Instituto Nacional 
do (Frio. 

9 — Disseminar as soluções de realização sim- 
ples a nível rural, usando materiais locais. 

10 — Lançar a secagem, a nível de exploração 
agrícola. Estabelecer a viabilidade econó- 
mica da dessalinização, da utilização de es- 
tufas, etc. 

11 — Lutar decididamente contra a separação das 
Escolas das realidades concretas da econo- 
mia do país. Apoiar decididamente quem 
se mostra capaz. Fornecer-lhes meios ade- 
quados, realisticamente, Lutar contra a pro- 
liferação de actividades descoordenadas, 
mas lutar também contra o centralismo as- 
fixiante e burocrático que estiola a criati- 
vidade e a inovação responsável, Tal cen- 
tralismo existiu e reforça-se pela via da 
concessão dos apoios financeiros e outros. 
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C — MÉDIO E LONGO PRAZO 


Poder-se-ia planificar uma actuação e estabelecer 
metas. Todavia, se o imediato e o curto prazo não 
forem vigorosamente encarados, de pouco vale falar 
do médio e longo prazo. Se o imediato e o curto prazo 
não forem vigorosamente encarados, no médio e longo 
prazo estaremos inteiramente colonizados e transfor- 
mados em fornecedores de mão-de-obra barata... e de 
capatazes diplomados. Os títeres não planeiam, em- 
bora possam votar. 

Em vez de falar no médio e longo prazo, é ful- 
cral assegurar que a actuação presente não vai liqui- 
dar o futuro. Não fazer nada é também um modo 
de actuação, 

A exposição anterior insere-se na linha de um 
futuro diferente, de um futuro de dignidade humana 
recuperada e vivida. O futuro começa hoje. É pois 
imprescindível que ataquemos as tarefas que ao Hoje 
cabem, em vez de nos dispersarmos na ilusão de um 
futuro que transfere para o milagre a solução de 
erros que não queremos enfrentar e corrigir, Hoje! 
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RESUMO 


Como os recursos para a investigação sempre 
escassos analisa-se a estratégia da sua repartição 
pelos sectores fundamentais a partir do seu efeito 
multiplicativo ou factor de impacto. Tal impacto tem 
em conta a estrutura actual da indústria mecânica 
do país e os factores essenciais que têm motivado 
o seu desenvolvimento. Finalmente mostra-se como 
a investigação em Termodinâmica Aplicada repre- 
senta um sector chave para um desenvolvimento 
moderno a caminho da progressiva automatização. 


1 — INTRODUÇÃO 


1.1. A primeira dificuldade em discutir o pla- 
neamento da investigação no sector da Engenharia 
Mecânica está em definir o que se entende por este 
sector. A Engenharia Mecânica é um sector essen- 
cialmente interdisciplinar, cobrindo um vastíssimo 
leque de aplicações. Entre nós, alguns dos seus domí- 
nios mais relevantes não entraram ainda, no senso 
comum como sendo do seu âmbito. Tal atitude 
decorre, essencialmente, do facto de à indústria 
mecânica se associar a manufactura, a produção, e 
se desconhecer a investigação básica que lhe é apli- 
cável e se traduz em projectos e realizações sem 
excessivo desenvolvimento intermédio. 

Na Engenharia Mecânica deve encontrar-se uma 
síntese dos seguintes sectores fundamentais: 


a) Termodinâmica Aplicada (Energia — sua 

transferência e conversão, sobretudo nos 

aspectos de calor e massa e mecânica dos 

fluidos). 

Identificação e cálculo das solicitações (for- 

cas) e das condições (temperatura, corrosão, 

etc.) em que o equipamento irá funcionar. 

c) Conhecimento e especificação dos materiais 
mais adequados à concretização do projecto. 

dy Tecnologia de fabricação propriamente dita. 


b 


e 
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SUMMARY 


Because research funds are always scarce, its 
allocation must be based on impact analysis. This 
analysis takes into account the present structure 
and development of the mechanical industries in 
Portugal and the main factor behind its development. 
It is show how research in Applied Thermody- 
namics is a key sector for modern development in 
a way to progressiva autonomy. 


1.2. Só a última é tradicionalmente considerada 
como especificamente mecânica, e dessa atitude se 
encontra reflexo nos programas dos próprios cursos 
de Engenharia Mecânica até há relativamente poucos 
anos. Nessa atitude se filia também o carácter arte- 
sanal de um grande sector, e a dependência tecnoló- 
gica em relação ao exterior do que nele constitui a 
fracção mais avançada. 

Ora, a indústria mecânica constitui uma infra- 
estrutura básica de qualquer país industrializado, e 
é factor essencial no seu desenvolvimento. No nosso 
país, em termos de postos de trabalho e de contri- 
buição para o P.N.B. a importância da indústria 
mecânica é muito apreciável e vai certamente cres- 
cer. Todavia, devemos acentuá-lo, o seu valor acres- 
centado está muito aquém do que poderia e deveria 
ser, porque a fracção que incorpora de matéria 
cinzenta é muito diminuta. Deste facto resulta que a 
sua capacidade de criação de emprego está muito 
longe do que poderia alcançar, e a valorização dada 
ao trabalho incorporado muito pequena. Do facto 
resultará que, embora globalmente o volume de 
emprego e de contribuição para o P.N.B. venha 
naturalmente a aumentar, esse aumento se verifique, 
sobretudo, por efeito de condições favoráveis de mão- 
-de-obra pouco qualificada e em consequência da 
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exploração dessa reserva por empresas estrangeiras 
ou suas subalternas nacionais ou nacionalizada. E se 
as condições do presente obrigam a que assim seja, 
a emancipação só poderá surgir se decididamente 
se valorizarem os factores a) b) c) d), pela ordem 
indicada, e simultaneamente se estimular e valorizar 
a sua intepenetração. 

As razões porque se advoga esta estratégia 
decorrem de uma análise do efeito multiplicativo do 
impacto de cada um deles. 


1.3. Deverá ainda acentuar-se que o financia- 
mento da investigação em Engenharia Mecânica é 
singularmente diminuto relativamente a outros ramos 
da Engenharia não apenas em valor absoluto como 
também, (e em situação muito mais agravada) 
quanto à percentagem do P.N.B. para que contribui. 
Fruto dos condicionalismos habituais, não pode deixar 
de se chamar veementemente a atenção para o facto 
de existirem três centros universitários neste domí- 
nio (dois em Lisboa e um no Porto) dois dos quais 
criados em pouco mais de um ano a quando da 
reestruturação dos Centros do hoje Instituto Nacio- 
nal de Investigação Científica (I.N.I.C.). 


O Centro de Termodinâmica Aplicada e Mecânica 
dos Fluidos, antigo Núcleo de Estudos de Engenharia 
Mecânica, cuja actividade se tem situado no domínio 
da investigação fundamental aplicável, e tem promo- 
vido a sua aplicação, apesar de mais de dez anos de 
existência, continua a não dispor de meios ou instala- 
ções que sejam sequer comparáveis aos que já hoje 
existem sobretudo no âmbito da Física e Química, 
como sejam o Complexo Interdisciplinar e o Instituto 
de Física e Matemática, ambos do I.N.I.C. 

Sem que se ponha em causa a necessidade de 
dotar aqueles domínios de meios adequados, nem por 
isso pode deixar de assinalar-se a necessidade de 
urgentemente serem corrigidas tais distorções de 
modo a estabelecer-se uma repartição de fundos 
mais consentânea com o impacto no desenvolvimento 
económico do país, e no seu potencial científico e 
tecnológico. 


2 — ANALISE DOS FACTORES DE IMPACTO 


2.1. A actual estrutura da indústria mecânica 
nacional é típica de um sector de fraco índice 
tecnológico cujo desenvolvimento se processou ao 
sabor de tradicionais factores de mercado, traduzido 
em milhares de pequenas unidades de implantação 
local que laboram como subsidiárias na manutenção 
de equipamento importado, ou no manufactura de 
produtos tradicionais (cutelaria, material agrícola, 
etc.), a par de um número muito restrito de grandes 
unidades inteiramente dependentes de tecnologia im- 
portada e de que são quase totalmente subsidiárias 
mesmo quanto à colocação de produtos. 


A expansão do sector irá naturalmente verifi- 
car-se com as características que possui, sendo espe- 
rável uma melhoria das tecnologia de fabrico pelo 
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simples efeito dos factores em jogo. Tal significa 
que não se prevê um grande impacto, ao nível das 
tecnologias de fabrico, de um significativo incre- 
mento da investigação aplicada nesse domínio, salvo 
em domínios muito restritos. Impacto importante, 
ao nível sobretudo da pequena e média empresa 
poderá sim verificar-se actuando ao nível da forma- 
ção profissional, dos cursos intensivos e da recicla- 
gem. Neste sector, é pois a actividade que possa 
estimular esta intervenção a que deve ser priorita- 
riamente contemplada. Tal não significa que a inves- 
tigação aplicada nesta zona seja descurada, mas 
sim que ela deve estimular basicamente uma activi- 
dade de formação. Nesta óptica o patrocínio de 
investigação nos domínios da soldadura e da fun- 
dição (ferro e metais não ferrosos, areias de fundi- 
ção, técnicas de moldagem, etc.), apresenta-se como 
a mais adequada. 


No escalão seguinte de impacto crescente situa-se 
o sector das máquinas-ferramentas, do seu controlo 
automático a vários níveis (peunático, hidráulico, 
digital, etc) e dos processos de enformação e corte 
de chapa. Deve, paralelamente, apoiar-se o desenvol- 
vimento da metrologia, não só pelo efeito geral na 
melhoria da qualidade (e uniformidade) da produção, 
como sobretudo pelo estímulo à mecânica de preci- 
são, que constitui um sector já importante e cujas 
potencialidades estão longe de se encontrar devida- 
mente exploradas. 


A filosofia que preside às recomendações ante- 
riores centra-se no facto de que a tecnologia do 
fabrico evolui com a solicitação e esta origina-se 
sobretudo no projecto e desenvolvimento de novos 
produtos e novos equipamentos. 


Nas linhas actualmente existentes no I.N.I.C€. 
é sobretudo o Centro da Universidade do Porto o que 
melhor se situa neste domínio. 

No concernente à soldadura refira-se o CEMUL, 
e as ligações existentes entre ele, o Instituto de 
Soldadura e o CIMAF. 


2.2. Num escalão qualitativo superior, situa-se 
o conhecimento e adequada especificação de mate- 
riais. Estreitamente ligado ao anterior, o carácter 
multiplicativo do investimento é maior desde que a 
disseminação dos resultados adquiridos seja assegu- 
rada. Este domínio tem estreita ligação com a meta- 
lurgia (tal como muitas técnicas de fabrico) e 
também com o vasto sector dos materiais não metá- 
licos; um e outro têm sido tradicionalmente tratados 
no âmbito da engenharia química, devendo por esse 
facto notar-se que a óptica de maior interesse para 
a Engenharia Mecânica é frequentemente descurada 
constituindo um bloqueio, raras vezes identificado 
na transferência para as aplicações de resultados 
obtidos pela Química ou Metalurgia. A óptica fun- 
damental para a Engenharia Mecânica é a das 
características mecânicas dos materiais, a sua ade- 
quação a finalidades específicas, a sua facilidade 
em ser trabalhada, a estabilidade dimensional, a 
durabilidade e envelhecimento, a da sua adaptação 
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através de tratamentos apropriados (tratamentos 
térmicos, por exemplo), etc. 
O CEMUL inclui, nas suas linhas de investi- 


gação, parte destes objectivos. 


2.3. A identificação e cálculo de solicitações 
situa-se na zona do impacto crescente. Não nos 
alongaremos sobre ele, por nos parecer demasiado 
óbvio. Acentuamos apenas que é imprescindível um 
decisivo apoio ao cálculo numérico, única via de 
exploração no grau mais elevado de novos conceitos 
de dimensionamento que optimizem simultaneamente 
a economia do material e a exploração máxima das 
suas características. Tal só será conseguido com o 
fácil acesso a meios poderosos de cálculo automático. 
A sua generalização, devidamente estruturada, é 
susceptível de efeito imediato, se forem simultanea- 
mente estimuladas duas frentes: o desenvolvimento 
de novos métodos de cálculo, seguros e eficazes; a 
normalização e extensa difusão de «programas de 
projecto», de tal modo concebidos que um utilizador 
corrente na indústria não seja obrigado a dominar 
os pormenores, permitindo a generalização de algo 
semelhante ao uso de um automóvel ou um telefone, 
em que é desnecessário conhecer a fundo como fun- 
cionam para eficazmente se poderem utilizar. 


No desenvolvimento de novos métodos poderão 
convergir várias linhas de investigação já existentes, 
devendo assinalar-se a importância do encontro útil 
de matemáticas e engenharias. Tanto neste desen- 
volvimento, como na sua disseminação, o reforço 
dos meios do Centro de Cálculo das Universidades 
de Lisboa, dadas as suas características de serviço, 
é fundamental, sobretudo no que se refere ao tele- 
processamento e ao tratamento óptico e gráfico da 
informação. 


O desevolvimento de tais métodos de cálculo 
centrou-se até agora, sobretudo, no Centro de Ter- 
modinâmica Aplicada, no C.C.U.L. e numa das 
linhas de investigação do CEMUL. 


3 — TERMODINÂMICA APLICADA 


3.1. A energia representa um sector vital em 
qualquer economia. Entre nós, tem existido a ati- 
tude simplista de a reduzir à hidroelectricidade e ao 
transporte e conversão de energia eléctrica. Com o 
aumento da componente térmica na produção de 
energia eléctrica, começou finalmente a emergir a 
ideia de que ele poderia ser bem mais vasto. Ora, 
se na hidroelectricidade e no transporte as perdas 
já são importantes, tal efeito é fortemente acen- 
tuado no caso das térmicas, sejam elas convencionais 
ou nucleares. Por outro lado, tanto na produção 
como na transformação da energia eléctrica, os 
aspectos referentes à transferência de calor são 
dominantes. Por sua vez, a transferência de calor e 
de matéria e a mecânica dos fluidos estão intima- 
mente associadas, e a unidade de todo o processo 
sobressai logo que é encarado um pouco acima do 
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empirismo tradicional. A sua vital importância e 
impacto tornam-se evidentes logo que o problema 
é colocado na perspectiva fundamental da termodi- 
nâmica. Como esta perspectiva vem sendo sistemati- 
camente adoptada pelo signatário, houve sempre o 
receio de a sobrevalorizar, acompanhada do senti- 
mento de para a sua importância se não ter devida- 
mente chamado a atenção, sobretudo num meio 
pouco receptivo e preparado para a compreender. 
Tal envolvimento pessoal deverá naturalmente ser 
tido em conta no que se segue, embora se procure 
nele o máximo de objectividade. 


3.2. Se examinarmos as indústrias mecânica e 
química de qualquer país desenvolvido, imediata- 
menet descortinamos que a transferência de calor e 
massa e a mecânica dos fluidos se encontra na raiz 
mesmo da maioria dos processos industriais na 
indústria química e na da energia, e que o seu efec- 
tivo domínio é um factor essencial para o efectivo 
domínio do projecto e da tecnologia de vastiíssimos 
sectores industriais. Na realização do seu equipa- 
mento ,a indústria mecânica tem papel fundamental. 
Na coordenação e concatenação deste equipamento 
diverso, desempenham importante papel as «Empre- 
sas de Engenharia». 


Pela complexidade intrínseca dos fenómenos 
físicos em causa, a parte empírica teve até há 
pouco papel dominante, e daí a sua compartimen- 
tação. O advento dos computadores e dos métodos 
de cálculo a que deu viabilidade permitiu um subs- 
tancial salto qualitativo no domínio unificado do 
seu tratamento e eumentou substancialmente a sua 
capacidade de impacto nos sectores mais variados. 
Esta foi a via seguida pelo CTAMFUL e a razão 
por que veio a desenvolver métodos numéricos e 
analíticos autónomos para a solução de sistemas de 
equações em derivadas parciais não lineares, e 
investigação fundamental de natureza experimental 
que lhe permitiram ou permitirão intervir em aspec- 
tos tão variados como o projecto de fornos; câmaras 
de combustão da indústria vidreira, cerâmica, dos 
cimentos, das centrais térmicas; da conservação de 
alimentos (indústria frigorífica), da hidrodinâmica 
naval, do aproveitamento da energia solar, geotér- 
mica, das ondas; no desenvolvimento dos métodos 
de previsão e dispersão de poluentes, etc. 

Por outro lado, este sector foi reconhecido pela 
UNESCO como dos de maior impacto no «take off» 
de países em vias de desenvolvimento na criação 
de uma tecnologia autónoma e capaz de extrair um 
máximo dos recursos naturais existentes nesses 
países. 

As razões anteriores ,a experiência já existente 
e a necessidade de identificar as zonas de impacto 
máximo num esforço concertado de investigação e 
desenvolvimento, sobretudo numa altura em que 
tantos mal entendidos parece existirem no modo 
como uma investigação fundamental de base pode 
contribuir, a curto prazo, para uma crescente eman- 
cipação tecnológica nacional, levam-nos a considerar 
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este sector como dos mais incisivamente capazes de 
ultrapassar múltiplas contradições. 

Porque não só o domínio da transferência de 
calor e massa e a mecânica dos fluidos são decisivos 
na realização de um laser químico, como no confi- 
namento e estabilidade de um plasma na fusão 
nuclear, como no projecto de um prosaico fogão de 
cozinha que aproveite ao máximo as calorias de um 
combustível cada vez mais raro e caro. Para o con- 
seguir, precisa incorporar os últimos avanços con- 
seguidos pelos matemáticos na demonstração de teo- 
remas de existência e unicidade, de equações em 
derivadas parciais em campos contínuos ou aleató- 
rios, como precisa da química física os valores de 
viscosidades ou condutibilidades térmicas de misturas 
binárias ou ternárias, de coeficientes de termo-difu- 
são, de velocidades de reacção; da física de películas 
finas as propriedades radiantes de superfícies; dos 
especialistas em electrónica a realização de conta- 
dores de hiper-frequências para a sua anemometria 
laser, etc. 

Consideramos pois que a este conjunto, generi- 
camente designado por Termodinâmica Aplicada, deve 
não apenas ser reconhecido o estatuto de Investi- 
gação Fundamental como o da Aplicada (que sabe 
como e onde aplicar) como lhe devem ser conce- 
didos os meios correspondentes a sector prioritário, 

Actualmente, neste sector, existe apenas um 
Centro em Lisboa, e uma linha de investigação no 
Centro de Engenharia Mecânica da Universidade do 
Porto. 

Para além de ser reconhecida prioridade à inves- 
tigação fundamental na transferência de calor e 
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massa e à mecânica dos fluidos, são de elevada 
importância as suas aplicações: 


— conservação de alimentos (indústria frigori- 
fica, de leofilização) 

— combustão industrial (caldeiras, fornos, etc.) 

— hidrodinâmica naval 

— aproveitamento de energia solar 

— aproveitamento de energia do vento 

— aproveitamento de energia das ondas 

— aproveitamento do calor perdido, no âmbito 
de um melhor aproveitamento da energia 

— motores de explosão e combustão interna 

— utilização e valorização de isolantes 

— exploração agro-pecuária 

— poluição 

— conforto da habitação. 


Em todos os sectores mencionados, existem já 
contactos ou projectos em curso com organismos 
públicos e privados, estando assegurada a utilização 
dos resultados que se obtiverem. Não existe, portanto, 
o problema da falta de interesse nem de aplicação 
de resultados. Existe sim o dos meios humanos e 
materiais nos Centros que o podem fazer. Cabe 
às entidades de quem depende a sua concretização, 
darem conteúdo à afirmação tantas vezes expendida 
de que a Universidade deve estar ao serviço do 
País. Quando, na Universidade houve centros que 
largamente ultrapassaram as realizações que seria 
legítimo esperar dos meios que lhe facultaram, é 
tempo de exigir que haja quem cumpra o que mini- 
mamente lhe compete. 
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Fenómenos de combustão com importância em Engenharia 


RESUMO 


Portugal depende dos combustiveis fósseis para 
satisfazer mais de 80% das suas necessidades ener- 
géticas presentes, ' 

Um conhecimento mais aprofundado da Mecá- 
nica dos Fluidos e da Química dos Sistemas com 
Combustão será indispensável se se pretender melho- 
rar o aproveitamento dos combustíveis e minorar os 
efeitos prejudiciais das emissões que lhe estão 
associados. 

A linha de Combustão do C.TAM.F.U.L. pro- 
põe-se contribuir para este esforço, desenvolvendo 
métodos de cálculo de câmaras de combustão e va- 
lidando-os através de experimentação laboratorial e 
de medidas em instalações industriais, divulgando os 
resultados obtidos e realizando cursos intensivos na 
sequência do que efectuou em Dezembro de 1975. 


INTRODUÇÃO 


Desejaria ajudar a situar o tema com a constata- 
ção de dois factos: 


1 — nos últimos 200 mil anos o Homem tem utiliza- 
do a combustão para os mais diversos fins e 


não teve à sua disposição nem modelos físicos. 


complexos nem computadores, excepto nos últi- 
mos 20 anos. 


2 -— as previsões optimistas indicam que 90% do 
combustível fóssil existente na Terra se esgo- 
tará dentro de 100 anos, (óleo e gás natural) 
e em 300 anos o carvão.(1) 


donde se pode inferir que os esforços adicionais no 
sentido de melhorar a compreensão actual dos fenó- 
menos de combustão serão irrelevantes. 

Em contrapartida faço notar os três seguintes 
aspectos: 


1 —o valor do combustível gasto anualmente no 
mundo é de cerca de 10” Escudos (p.e. Portugal 
consumiu em 1974 mais de 10” Escudos) ou se- 
ja, a humanidade dispende em energia uma 
parte substancial das suas possibilidades, 


MÁRIO NERY R. NINA (*) 


ABSTRACT 


Over 80% of the energy used in Portugal today 
is derived from fossil fuels. 

A better understanding of the Fluid Mechanics 
and the Chemistry of Reacting Systems will be re- 
quired in order to improve the efficiency of combus- 
tion processes and simultaneously decrease the 
amount of harmful effluents. 

The Combustion Group of CTAM.F.U.L. is 
concentrating efforts in the development of better 
furnace prediction methods, validated through ezx- 
perimental data obtained both in the laboratory and 
in industrial instalations, and offering short specia- 
lised courses like the one that took place in Decem- 
ber 1975. 


2— À agressão ao ambiente provocada pela utiliza- 
ção de tão elevada quantidade de combustível 
fóssil tem consequências diversas e pesadas. 


3 — tanto quanto se pode prever, as fontes alterna- 
tivas de energia (solar, geotérmica, das marés, 
etc.) apresentam problemas de armazenazem 
e/ou transporte. Por exemplo, no caso específico 
dos transportes aéreos e rodoviários não está 
à vista a viabilidade de substituição da energia 
química dos combustíveis líquidos fósseis ou não 
por outra forma alternativa. 


A procura de soluções para os problemas que 
decorrem directamente dos aspectos atrás enuncia- 
dos deverá orientar-se, essencialmente, em duas vias 
fundamentais e distintas: 


— definição de uma política energética global (na- 
cional) no sentido de uma racionalização produ- 
ção-utilização, redução dos desperdícios, das emis- 
sões, controlo das agressões ao meio ambiente, 
com apoio de legislação e meios de intervenção 
adequados; é 


(') Estimativas dependentes da tecnologia e dos respectivos custos. 
(*) Assistente do 1. S..T., Bolseiro do CT. A. M. F. U. L. — Centro de Termodinâmica Aplicada e Mecânica de Fluídos 


das Universidades de Lisboa. 
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— progresso científico e tecnológico de modo a me- 
lhorar a relação rendimento-custo dos sistemas e 
controlar a poluição gerada pela combustão e isto 
tendo em conta que a qualidade dos combustíveis 
poderá piorar, pelo menos, a médio prazo. 


Limitando-me a este segundo aspecto de melho- 
ria tecnológica e retomando o que atrás disse ao 
tentar dar pontos de apoio para localizar o tema, 
poderei afirmar que se reveste de particular impor- 
tância a procura continuada de uma melhor com- 
preensão dos fenómenos de combustão porque: 


— um pequeno aumento de rendimento conduz a 
economia de verbas substanciais; 


— a poluição tem necessariamente de ser controlada 
dado que os custos da agressão que provoca são 
enormes, a qualidade dos combustíveis continuará 
a piorar e a quantidade utilizada e o preço deve- 
rão aumentar pelo menos no futuro próximo; 


— os combustíveis líquidos serão insubstituíveis du- 
rante vários séculos em sectores importantes da 
actividade humana e a tecnologia actual de com- 
bustão de hidrogénio (combustível do futuro) está 
longe de corresponder às exigências desse futuro 
(p.e. altos níveis de emissão de NO», limites de 
detonação inadequados, baixa eficiência dos sis- 
temas de controlo de diluição nas câmaras de 
combustão). 


EXEMPLOS DE COMBUSTÃO NA INDUSTRIA 


Passarei agora em revista exemplos de fenóme- 
nos de combustão tentando localizar as áreas onde 
a investigação em combustão tem interesse para o 
Engenheiro, fabricante de equipamento e utilizador. 


Fig. 1-—- Forno-tanque de fusão de vidro (corte) 


A — o forno industrial 


A figura 1 mostra um corte esquemático de um 
forno-tanque de fusão de vidro. A diversidade de 
tipos de fornos utilizados na indústria é grande tor- 


188 


nando este equipamento especialmente interessante 
do ponto de vista de emprego de métodos de cálculo. 
Os custos de capital são elevados o que obriga a 
uma grande prudência na inovação ou até na simples 
extrapolação. Por exemplo, conseguir uma distribui- 
ção adequada do fluxo de calor por radiação é im- 
portante pois permite baixar os custos de investi- 
mento (permitindo o emprego de materiais menos 
dispendiosos) e aumentar o rendimento térmico e a 
produção unitária, 

A conversão de equipamento existente para ou- 
tro combustível (passagem de carvão a fuel-oil ou 
gás natural) pode ser analisada com o auxílio dos 
métodos de simulação já existentes. 

Resumindo, enunciam-se algumas das exigências 
correntes: 


— distribuição de temperatura necessária ao proces- 
so, por conveniente repartição entre transmissão 
por radiação e convecção; 


— combustão completa com emissões dentro dos ní- 
veis eventualmente impostos por lei; 


— estabilidade de chama e ausência de oscilações 
e ruído; 


— tolerância às flutuações das características do 
combustivel disponível no mercado; 


— rendimento térmico global adequado ao preço 
previsível do combustível para a vida da instala- 
ção e ao capital a investir. 


B-— a câmara de combustão dos geradores de vapor 
(figura 2) 


Este órgão fundamental na produção de energia 
tem sido objecto de cuidadosos e variados estudos 
e tem de tal modo atraído as atenções dos engenhei- 
ros de combustão que dominou a cena durante anos. 
Os problemas que kevanta são muito importantes e 
complexos mas não devemos esquecer que represen- 
ta apenas um tipo particular de câmara de combus- 
tão que classificarei como combustão de alta inten- 
sidade (pois uma das formas de reduzir o custo ini- 
cial é diminuir o volume da câmara). 

Em oposição a esta situação está o exemplo 
anterior do forno-tanque de fusão de vidro que tem 
«chamas longas», No aspecto da estabilidade de cha- 
ma e da necessidade de equipamento de detecção 
este problema é peculiar às chamas de alta intensi- 
dade em câmara de paredes frias. 

Por outro lado, as dimensões das fornalhas dos 
grandes geradores de vapor para Centrais Térmicas 
(600 — 1000 Mw) e a tendência para reduzir o núme- 
ro de queimadores cria condições faviráveis à exis- 
tência de fenómenos de deslocação de chama que 
provocam assimetrias nos fluxos de calor com graves 
consequências tanto na operação como na vida do 
equipamento. 
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Fig. 2—- Câmara de combustão de gerador de vapor (600 Mw) 


Algumas características desejáveis são em re- 
sumo: 


— combustão completa com níveis de emissão con- 
trolados; 


— controlo das corrosões que afectam os órgãos sub- 
metidos a temperatura elevada; 


— ausência de ruído e oscilações; 


— bom comportamento com a variação de carga em 
especial temperatura de vapor constante; 


— boa estabilidade de chama (chama de alta inten- 
sidade em câmara de paredes frias). 
C—- o motor de combustão interna 


Os problemas que se nos deparam neste tipo de 
sistema são variados e incluem a procura de maio- 
res potências e rendimentos e, mais recentemente, 
a limitação das emissões consideradas nocivas. 
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No caso do motor diesel aproveita-se a ignição 
espontânea da mistura inicial sendo depois a com- 
bustão controlada pelo ritmo da injecção do combus- 
tível. A busca de mais potência conduz a usar todo 
o oxigénio disponível na câmara e como consequên- 
cia à emissão de fumo; procura-se ainda tornar si- 
multaneamente a combustão tão adiabática e silen- 
ciosa quanto possível, 

O motor a gasolina tem o seu rendimento essen- 
cialmente limitado pelo fenómeno de auto-ignição 
que é a todos os títulos prejudicial e que é contro- 
lado pelos construtores através de geometrias ade- 
quadas para a câmara e por meio de aditivos no 
combustível, mas o problema continua em aberto 
pois tem interesse conseguir melhores rendimentos 
usando combustíveis de mais baixo índice de octano. 

Como características e problemas típicos do 
diesel temos: 


— vaporização e/ou combustão de gotas; 


— velocidade limitada pela dinâmica das sucessivas 
fases do processo de combustão; 


— emissão de fumo e ruído; 
- arranque a frio difícil. 
No caso do motor a gasolina: 


— à cinética das reacções químicas condiciona o 
avanço da frente de chama no sistema com pré- 
-mistura;, seria desejável conseguir câmaras de 
combustão funcionando com misturas pobres; 


— controlar as emissões de hidrocarbonetos, monó- 
xido de carbono e de óxidos de azoto. 


D—a câmara de combustão de turbina de gás de 
aviação (figura 3) 


Pretende-se um sistema de queima compacto, 
oferecendo baixa resistência à passagem dos gases 
(devido ao efeito sobre o rendimento do ciclo) e 
fornecendo à turbina uma corrente de gás a tempe- 
ratura espacialmente uniforme. 

Do ponto de vista de emissões, o problema pre- 
sente não é grave pois existem boas condições de 
dispersão excepto na fase de manobra nos aeropor- 
tos. Assim a optimização do sistema de combustão 
está condicionada por vários parâmetros irreconciliá- 
veis e aqui também os métodos de cálculo poderão 
vir a substituir os ensaios em protótipo: 


— combustão completa e estável; 
— baixa perda de carga; 


— temperatura uniforme à saída (tensões térmicas 
nas pás da turbina). 
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UOMPRESSOR 


Fig. 3 —- Câmara de combustão de turbma de gds 
(zona de chama; 


Estes quatro exemplos servem para ilustrar apli- 
cações industriais onde um melhor conhecimento 
dos fenómenos de combustão pode traduzir-se em 
benefícios de natureza diversa, tais como, maior 
economia, maior potência, menor ruído e poluição, 
vida mais longa para os componentes do sistema, 
que são outras tantas preocupações do projectista 
e do utilizador. 


PERSPECTIVAS DA INVESTIGAÇÃO EM COMBUSTÃO 


Portugal recorre aos combustíveis fósseis para 
satisfazer cerca de 80% das suas necessidades ener- 
géticas, situação que não virá a ser substancialmen- 
te alterada a médio prazo. Por outro lado, o país 
importa, presentemente, a maioria dos equipamentos 
onde é levada a efeito a queima do combustível 
(quer sob a forma de aquisição do projecto quer por 
directa importação de equipamento) no entanto, a 
área de intervenção restrita que está à disposição 
dos técnicos portugueses é importante, quer na defi- 
nição das características e escolha dos equipamentos 
quer na sua exploração e melhoria, e não tem sido 
mesmo assim totalmente utilizada. 

O desenvolvimento e permanente actualização 
de um corpo de conhecimentos capaz de permitir 
atacar eficazmente os problemas dos sistemas de 
combustão com que se defrontam os técnicos portu- 
sueses terá diversas repercussões e a curto prazo 
dará contribuição positiva para o problema mais 
geral da Energia e da sua racional utilização. 

Sem pretender apresentar exaustivamente tópi- 
cos de combustão onde existem áreas por clarificar, 
enunciam-se alguns temas que são objecto de aten- 
ção em programas de trabalho no C.'T.A.M.F.U.L.: 


1 —escoamentos turbulentos com reacção quimica — 
as grandes chamas industriais estão aqui incluí- 
das e existem modelos para descrever a transfe- 
rência de energia, de massa e os aspectos glo- 
bais da cinética química que intervêm, O aper- 
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feiçoamento destes modelos depende criticamen- 
te de experimentação adequada. 


2 — instrumentação —o projecto e construção de 
aparelhagem capaz de permitir clarificar e ava- 
liar características particulares dos escoamen- 
tos com reacção é essencial; em particular, o 
desenvolvimento do anemómetro laser possibili- 
tará a observação dos mecanismos de turbulên- 
cia relevantes, 


3 — limites de ignição em misturas ricas e po- 
bres e limites de extinção — as misturas afasta- 
das da composição estaquiométrica são mal co- 
nhecidas embora se revistam de grande impor- 
tância nas aplicações. 


4 — combustão de «sprays» — na grande maioria 
dos casos o combustível líquido é previamente 
pulverizado e os mecanismos que controlam a 
inter-acção das gotas com a atmosfera oxidante 
envolvente são mal conhecidos. O nível de emis- 
são de poluentes e a eficiência da câmara de 
combustão estão ligados aqueles mecanismos. 


5 — formação e destruição de particulas sólidas — 
é impossível calcular com rigor a transferência 
de calor por radiação numa chama luminosa 
desconhecendo os parâmetros que afectam a 
formação e o desaparecimento das partículas só- 
lidas. Estas partículas, por outro lado, vão con- 
tribuir para as emissões e são responsáveis pelas 
corrosões que se dão em áreas localizadas dos 
equipamentos. E pois importante vir a dominar 
estes processos que ocorrem em muitos sistemas 
de combustão com particular relevo para os que 
queimam fuel-oil e carvão. 


6 — métodos de cálculo para câmaras de combustão — 
todos os aspectos anteriores contribuem para o 
sucesso ou insucesso dos métodos de cálculo 
com os quais se pretende simular as condições 
existentes numa câmara de combustão com uma 
geometria e um combustível impostos. 


Estes métodos estão operacionais com as limita- 
ções que decorrem quer do incompleto conhecimento 
de certos mecanismos da mecânica dos fluidos e da 
cinética das reacções químicas quer da capacidade 
dos meios de cálculo à disposição do Centro, 


Dentro dos objectivos do CT.A.M.FU.L. e em 
particular da sua linha de Combustão inclui-se a 
divulgação do trabalho produzido, através de relató- 
rios de actividades periódicos, da publicação de arti- 
gos em revistas da especialidade e ainda do trabalho 
docente dos seus colaboradores no âmbito das Uni- 
versidades de Lisboa e mais especificamente em 
cursos intensivos abertos à participação da indústria, 
como o que se realizou en Dezembro de 1975. 
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A calorimetric experimental furnace 1) 


RESUMO 


Projectou-se e construiu-se uma fornalha de sec- 
ção quadrada (500X500 mm) com o objectivo de per- 
mitir obter valores experimentais para validação 
de métodos de cálculo de câmaras de combustão 
e fazer ensquios de queimadores e de corrosão de alta 
temperatura, 

4 estrutura da fornalha permite variações na 
geometria dos circuitos dos gases e da úgua de 
arrefecimento. 

Apresentam-se resultados de ensquios de um quei- 
mador do tipo cone rotativo, onde se mostra a in- 
fluência da carga e do excesso de ar na distribuição 
do fluxo de calor. 


INTRODUCTION 


Although the general laws of Combustion and 
Heat Transfer are well understood, fuel-fired furnace 
design remains an art in many respects; the need to 
reduce costs and to control emission of harmful 
combustion products are two important reasons to 
look for better and more accurate design methods. 

Within the scope of the joint research program 
between Division of Applied Thermodynamices of 
I S.'T. and L. E. M, it was decided to design and 
built an experimental calorimetric furnace to be 
used primarily for testing of prediction procedures 
applied to furnace design and more specifically the 
prediction of temperature distributions, and to help in 
problems directly connected with the needs of local 
industry in furnace design, manufacturing and ope- 
ration. 


EQUIPMENT 


The design of the furnace body structure enables 
the easy modification of its geometry. The basic 
constitutive element, fig. 1, is a square ring, 500 mm 
by 500 mm, 100 mm thick. 


3.13. D. DOMINGOS, 
E. MANSO () 
M. N. R. NINA 


ABSTRACT 


A 500% 500 milimeter square section furnace was 
designed and built to provide data for the testing of 
prediction procedures apphed to jfurnace design and 
to enable experiments to be carried on high tempera- 
ture corrosion and burner design. 

The furmnace structure allows the geometry of 
the gas passage and the cooling circuit to be assem- 
bled in several ways. 

Some preliminary test results are presented and 
show the effect of excess of air and load on local 
heat flux using a rotating cup type of burner. 


A number of these elements can be juxtaposed 
in a way that there is continuity both on the gas and 
the water passages, Provisions were made to allow 
the water circuit to be obtained either by connecting 
the elements in series, in parallel-series or as sepa- 
rate groups of elements, 

Each element is self contained in what regards 
thermal insulation. A thermocouple well in the water 
circuit and two flame observation windows are pro- 
vided in each standard element. There are also non- 
-cooled refractory elements for specific uses, p.e. with 
gas recirculation inlet ports, with instrument doors. 

A. 90º refractory corner makes it possible to have 
a water cooled vertical section or a direct connection 
to the stack. 

The furnace is instrumented to allow measure- 
ments of metal temperatures, water temperatures, 
pressures and flow rates, gas temperatures (with an 
aspirating pyrometer) and combustion gases compo- 
sition, The furnace cooling water runs in closed cir- 
cuit passing through two sets of tubular heat 
exchangers. Two independent water flowmeters are 
installed, besides a measurement by weigth. The 
combustion air is supplied to the burner by an elec- 
tric motor driven centrifugal fan; two separate ori- 
fice meters (in accordance with B. S. 1042) allow 
the measurement of primary and secondary air flows. 


(1) Paper for the European Symposium on Combustion, Sheffield, September 1973. 


(2) LEM — Direcção-Geral de Combustíveis, Lisbon. 
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Figure 3 -— Overall view of the furnace 
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Figure £ -—- Cooling circuit within a single element 
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The fuel is measured by a weighting method. 
Its temperature is controled in two steps to keep a 
constant value at the bumner outlet. 

The cooling water temperature is measured at 
the inlet and outlet of each element through ther- 
mocouples inserted in a well around which the flow 
section is substantially reduced to minimize thermal 
lag and conduction errors. The thermocouples are 
connected to one automatic potentiometer and from 
time to time checked against a precision instrument. 
The thermocouples are connected in a differential way 
and the global error is checked against the overall 
temperature difference measured between the inlet 
and outlet temperatures of the cooling water. Gas 
temperature and composition are measured at the 
outlet section after a mixing bafle. A continuous 
sampling and recording system is used for CO, 
CO, O,, SO. 


7 


O 


1 10 20 nº of afomest JO 
E DDT RS 5 TT 25 177,08 
(8) 1000 2000 3000 mm 


Figure 4 — Influence of Load on local heat flux 


PRELIMINARY TESTS AND RESEARCH PROGRAM 


Trial tests have been conducted during several 
months to calibrate the instruments, improve the 
heat balance and to get familiarity with the dynamics 
of the furnace. The total heat balance according to 
the ASME boiler test code is usually closed within 
4-5 % which accounts for radiation and unaccount- 
ed losses. Most of the tests have been conducted 
with an automatic cup atomizing burner of the 
type commercially available for tube fired boilers. 
The aim was to measure the heat flux distribution 
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versus excess air and load. Fig. 4, 5 show some re- 
sults. These are plots in a conventional way, showing 
in the abcissa the length along the furnace and in 
ordinate a non dimensional heat flux expressed as 
a percentage of total heat imput, This type of burner 
showed little flexibility in what concerns load and 
low excess air. Im this test the whole surface is 
cooled and the refractory cone which should surround 
the burner according to its manufacturer was not 
used. Due to the bad atomizing characteristics and 
unfavorable flow field some unburned fuel accumula- 
ted in the bottom of the furnace which is reflected 
in the error of closure of the heat balance. 

The effect of swirl on pressure jet flames is 
presently under study and will be reported later. 

The experimental results got in this furnace 
will be used in connection with a basic research pro- 
ject in the turbulence structure of swirling jets and 
computer programs developed for boundary layers, 
and two and three dimensional flows coupled with 
heat and mass transfer and chemical reaction 
(Ref. 1, 4). In a separate paper it is shown how 
droplet trajectories will be predicted and used 
(Ref. 2). 

Due to the flexibility in design of the furnace, 
corrosion studies and heat transfer to convection 
banks are easily dealt with. Differences in tempera- 
ture between sections and its influence on heat 
transfer by radiation are feasible because each sec- 
tion can be independently cooled. In connection with 
this application it is probably worth mentioning the 
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Figure 5 — Influence of excess of air on local heat flux 
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problems which can arise with heated circuits either 
in series or parallel if the main pressure drop is at 
the outlet or in the heated parts themselves: for 
certain flow regimes the flow becomes unstable. 'The 
phenomena is well explained even in a pure static 
way (Ref. 3); when considered dynamically, new 
modes of instability occur. This furnace can show 
the phenomena easily in a controled way varying 
the inlet and outlet localized pressure drops, even 
without vaporization in the circuit. This type of 
phenomena has been considered at the design stage 
both to avoid it in heat transfer studies and also 
to promote it when necessary to investigate dyna- 
mic characteristics of the cooling flow, 
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Transferência de calor em regime estacionário em chaminés 
(Análise e Resolução Numérica) 


RESUMO 


O presente trabalho surgiu da necessidade do 
conhecimento da temperatura e velocidade dos gases 
à saída duma chaminé conhecidos os seus valores 
na base. Tais valores são importantes para a pre- 
visão da dispersão de poluentes emitidos pelas cha- 
minés para a atmosfera. Analisa-se a transferência 
de calor entre os gases e o exterior das chaminés. 
Calculam-se os coeficientes de transmissão super- 
ficial de calor e as propriedades do ar a 760 mmHg. 
Realizou-se assim um programa para cálculo da 
temperatura dos gases ao longo da chaminé, da sua 
velocidade na base e no topo e os perfis de tempe- 
ratura nas paredes interior e exterior. Analisam-se 
seguidamente alguns casos reais, sendo tratado o 
problema da corrosão de chaminés metálicas. 


1. INTRODUÇÃO 


Nos estudos de dispersão de poluentes emitidos 
por chaminés para a atmosfera dois dos dados 
essenciais numa situação prática concreta são o 
conhecimento da temperatura e da velocidade dos 
gases à saída da chaminé. A medição directa do 
valor da temperatura é em muitos casos difícil ou 
mesmo impraticável, e raras são as instalações que 
se encontram preparadas para o fazer. Tanto quanto 
é do nosso conhecimento não existe neste momento 
no nosso país nenhuma instalação com essa possi- 
bilidade de medida. 

Por outro lado a temperatura de saída dos gases 
da chaminé é dependente de vários factores atmosfé- 
ricos como a temperatura ambiente e a velocidade 
do vento. Estes factores não são constantes embora 
deles se possam tomar valores médios. Um outro 
factor que afecta grandemente aquela temperatura 
é o caudal dos gases expelido pela chaminé. 
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ABSTRACT 


The prevision of the dispersion of airborne 
chimney or stack effluents requires that the tempe- 
rature and the speed of the gases at the top of the 
stack should be known. This paper presents a way 
of computing those values given the temperature at 
the bottom, the mass flow rate of the gases and the 
caracteristic stack dimensions. An analysis of the 
heat transfer between the flue gases and the outside 
of the stack is made. The heat transfer coeficients 
and the properties of air at 29.92 inches of Hg are 
calculated. A computer program was developped. 
This program computes the longitudinal tempera- 
ture profiles of gases and of the inside and the 
outside stack surfaces, and the gas speed at the 
top and bottom of stacks. The last part of this paper 
analyses a few real cases and presents an example 
of prevision of minimum inside stack surface tempe- 
rature, an important parameter for preventing acid 
corrosion in metalic stacks. 


Os ábacos que existem para o cálculo da tempe- 
ratura de saída são construídos em função da altura 
e diâmetro da chaminé e da temperatura a que os 
gases nela entram. Qualquer variação da espessura 
das paredes altera as condições de transmissão de 
calor e portanto a temperatura que nos interessa. 

Para a resolução do problema foi desenvolvido 
um programa em linguagem BASIC para o calcu- 
lador Hewlett-Packard 98304. O programa pode ser 
obtido no Centro de Termodinâmica Aplicada e Me- 
cânica dos Fluidos das Universidades de Lisboa 
(antigo Núcleo de Estudos de Engenharia Mecânica). 
Este programa permitiu obter resultados respeitantes 
a uma chaminé de betão que se compararam com 
outras fontes. Foi também possível abordar o pro- 
blema da determinação da temperatura mínima da 
parede interna duma chaminé metálica. 


195 


2. CONSIDERAÇÕES GERAIS 


2.1. Sobre Transmissão de Calor 


O problema da transmissão de calor através e 
ao longo das paredes duma chaminé é fundamental- 
mente um problema bidimensional visto que supore- 
mos que existe uma simetria radial. Esta hipótese 
é feita por se desprezarem os efeitos da insolação 
e se supor que a atmosfera nas imediações da cha- 
miné se encontra toda à mesma temperatura. 

O regime de transmissão de calor é suposto 
permanentemente e portanto, teremos para a cha- 
miné 

GEE A mm 1 oºT qT 


dr? r or r2 090º dZ2 


=0 
ST 
que devido à condição de simetria radial Rd é 
O 


0 


92T 41 aT PRE 
dr2 r or dg 


Se o gradiente axial de temperatura for aproxima- 


2 
damente constante, o termo Ca pode ser despre- 


zável em relação aos outros termos. Pelos resultados 
obtidos podemos concluir à posteriori que a hispótese 
é válida. Assim, em geral, é 


o2T + 1 07. 
dr? r dr 


pelo que o problema pode ser considerado unidimen- 
sional. 

Como condições fronteiras teremos no interior o 
fluído constituído pelos gases a uma determinada 
temperatura e no exterior o ar atmosférico a outra 
temperatura. 


2.2. Natureza dos gases 


Estequiometricamente uma combustão exige 
quantidades fixas de combustível e comburente, que 
no caso que estamos a tratar é o oxigénio. Na prá- 
tica, no entanto, não se fornece ao combustível só 
a quantidade de ar que contém exactamente a quan- 
tidade de oxigénio necessária à combustão com- 
pleta duma determinada quantidade de combustível, 
pois para concentrações de oxigénio inferiores a um 
dado valor é muito difícil fazê-lo chegar à zona onde 
se está dando a combustão. Assim, torna-se impe- 
rativo que exista um excesso de ar, o que acarreta 
a consequência de nem todo o oxigénio ser substi- 
tuído por CO, e pelos outros produtos menores da 
combustão. A concentração destes produtos no ar 
diminui pois, com o aumento de excesso de ar. 

Tomando como exemplo um teor de CO, de 10 % 
podemos comparar o que acontece à composição 
gravimétrica, 
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Composição normal Com 10 % CO; 


Peso Fracção Peso Fracção Peso 
molecular molar molecular molar molecular 

Parcial parcial 

O, 32,000 0.2095 6,704 0,1098 3,914 
N, 28,016 0,7809 21,878 0,7809 21,878 
A 39,944 0,0093 0,371 0,0093 0,371 
Co, 44,01 0,0003 0,013 0,1000 4,401 
1,0000 28,966 1,000 30,164 


Como se poderá concluir os gases que neste caso 
se escoariam pela chaminé pouco diferem do ar 
quanto a propriedades físicas globais pelo que, com 
boa aproximação, poderemos considerar o efluente 
com as propriedades físicas do ar. Esta hipótese 
será adoptada em todo o trabalho. 


3. ANALISE DA TRANSMISSÃO DE CALOR 
3.1. Algumas notas preliminares 


Esta secção trata o problema de transmissão 
de calor entre os gases quentes que circulam no 
interior da chaminé e a atmosfera exterior a esta. 
Esta análise é efectuada por meio de balanços tér- 
micos nos subsistemas locais do sistema total que 
é a chaminé. Este processo foi adoptado devido ao 
escoamento interior da chaminé ser turbulento o que 
torna impossível um tratamento fácil das equações 
de movimento e consequentemente a determinação 
de perfis de temperatura. Para evitar tais inconve- 
nientes convencionou-se adoptar para temperatura 
dos gases a sua temperatura de mistura adiabática 
(bulk temperature) T, que em escoamentos interio- 
res a condutas de secção circular e supondo conhe- 
cidos os perfis de velocidade e temperatura é, 


“dx PR 
EV “Fs | | v 4riTirirdrd 
o pad h. o: U . O A 
b O v a “Or *R 
7 | | v (rirdrd) 
() 0) 
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As determinações dos coeficientes de transmis- 
são de calor e das propriedades dos gases são supos- 
tas completamente definidas pelo que o seu trata- 
mento é deixado para secções posteriores. 

Os regimes do escoamento e de transmissão de 
calor são considerados permanentes. 


3.2. Temperaturas de referência 


Como as propriedades dos gases são fortemente 
dependentes da temperatura e esta varia de ponto 
para ponto é necessário definir a temperatura a que 
aquelas se deverão determinar. 

Quando se trate de propriedades envolvidas na 
transferência de calor de gases para um meio sólido 
tomaremos como temperatura de referência a média 
das temperaturas que caracterizam o processo ou a 
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temperatura de mistura. Assim entre os gases da 
chaminé e a parede interior teremos 


Ali f — E 
sendo 


Ti,., — temperatura de referência interior 
T, — temperatura de homogeneização dos ga- 
ses da chaminé 


Para a transferência de calor entre a parede 
exterior e a atmosfera a temperatura de referência 
local é 

E PL 
oo e 
TO,er — 
2 
sendo 


Te,., — temperatura do filme exterior 
To — temperatura da atmosfera não afectada 
Te —temperatura da parede exterior 


Quando for necessário calcular propriedades do 
fluído características do escoamento no interior da 
chaminé a temperatura a utilizar será também a 


temperatura de homegeneização (T, ). 


3.3. Significado das grandezas médias 


Este parágrafo tem cabimento uma vez que o 
escoamento interior não é um escoamento isotérmico 
através de secções constantes, pois em chaminés a 
forma é geralmente cónica. 

Uma vez que se trata dum escoamento perma- 
nente o caudal mássico de gases é constante em 
qualquer secção que se considere. Vamos considerar 
duas quaisquer delas que sejam suficientemente pró- 
ximas (ver fig. 1). 


6 


| 
| 


figura 1 


Uma unidade elementar de massa de gases dM 
percorre o trajecto entre as secções 1 e 2, tendo 
nestas secções as propriedades que na figura se 
indicam. 
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Supondo que a variação de qualquer tempera- 
tura (a da parede interior, por exemplo) é linear 
ao longo de cada secção é imediato que os valores 
médios se verifiquem na secção média. Portanto as 
grandezas que vão entrar no cálculo dos balanços 
térmicos são grandezas que se referem à secção 
média de cada um dos troços. 


3.4. Balanço térmico 


Tomando uma secção elementar duma chaminé, 
obtida pela intersecção desta com dois planos para- 
lelos a uma distância elementar um do outro e 
perpendiculares ao eixo da chaminé (veja-se a fig. 1) 
obtemos o que se pode considerar um tubo cilíndrico 
oco de raio interior r, e raio exterior r,. 

O tubo é composto por n camadas concêntricas 
de condutibilidades térmicas K,, K,,..., K,. Os raios 
intermédios são definidos do modo como a figura 2 
mostra. Desprezaremos as resistências de contacto. 


(2) ln) 
RA 
K, Kn 
Fien4) 
| h = 
figura 2 


Sendo T, e T, respectivamente a temperatura 
da superfície interior e a temperatura da superfície 
exterior e com T, maior que T,, o fluxo de calor 
que passa do interior para o exterior numa secção 
de altura AZ é 


27 (TT, AZ 
q= peter 
io o q Ejo 1 Ea 
—+ + +... + In —— 
K, Fr, K, F, K, F, tn-1) 
(3.1) 


Dos gases que circulam na chaminé para a pa- 
rede interior o fluxo de calor é devido à acção 
conjunta da convecção e da radiação. Supondo que 
os gases se encontram à temperatura média T,, atrás 
definida, a expressão do fluxo é da forma 


q=2rr, 42 [b(T,—T)+e «(TT (3.2) 
sendo 


h, — coeficiente de convecção entre os gases 
e a parede interior 
e, — emissividade da superfície interior 
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O fluxo de calor entre a parede exterior e a 
atmosfera deve-se também à convecção e à radiação 
e a sua expressão é 


q=2rr, AZ [bh (T—T o) +e, (T—TID)] (33) 


sendo 
T,. — temperatura do ar atmosférico não per- 
turbado 
h, — coeficiente de convecção entre a parede 
exterior e a atmosfera 
e. — emissividade da superfície exterior 


A determinação dos coeficientes h, e h, é tra- 
tada na secção 4. 

Ao atravessar o troço de tubo os gases sofrem 
um abaixamento de temperatura AT. Denotando por 
T,, à temperatura de referência à entrada do troço 
elementar e continuando a supor que a variação de 
temperatura ao longo do troço é linear, a tempera- 


tura média de referéncia ER é 


a (3.4) 
2 


Significado análogo a A se atribui a T e T,. 
Ao atravessar então o troço os gases cedem a 
quantidade de calor 


1 
2 
= par AZ CraT Xx-—— 
q p 1 b At 


p é Cr são determinados à temperatura média de 


referência To e àt o tempo necessário para renovar 


completamente os gases que se encontram no interior 
do troço elementar considerado. Assim é 


AZ 
dt = -— 
V 


sendo V a velocidade média dos gases no troço. 
A equação toma a forma 


V 
- 2, aii a 
q — ade ada on TS cd de 1 a Rod 


por seu lado o caudal mássico dos gases é 


G = pVr E 
ou seja 

V= ds 

pT Fr, 

donde 

q = C,GAT 

(3.5) 
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Os valores de q em (3.1), (3.2), (3.3) e (3.0) são 
iguais pelo que podemos formar o seguinte sistema 
de equações 


4 mm 
277,82 [h,(T,— To) te (T— Ti)]= 


| = E E 
[277,42 [h, (1, —T) + eo (T;— T9]= 


o SUAR E) 
A 
27(T, E ih 
C GAT, = AZ 
A 
1 Xi 1 Fo 1 r 
sendo Ai=——lIn + —— In +...+t—— ln 
1 4 o H Ty n Fina) 
(3.6) 
ou seja 
T—T 
ES —4 4 1 2 
h(T,— T) + e19 na = 0 
T,—T 
4 4 1 2 
( h, E 8) + E99 O a ii 0 
2 (1, — To) 
e G GD E = 0 


Introduzindo agora a equação (3.4) ter-se-á 


an sT aT 
b b 4 4 
==" —EPao [5 — E] — 
2(T, RS 
——w———— = 0 
( A 2(T, = a) 
h(T—T, D+ ego (TT, d—————=0 
2n(T,— T,) 
a SEA, ee re A O) 


(3.7) 


que é um sistema algébrico não linear de 3 equações 
a 3 incógnitas (T,, T, e áT,). 


Neste sistema existem várias funções que depen- 
dem das temperaturas, que são: 
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at, 


h =h,(T9)=h (T—— 


) 


ST 


cume b 
C=C(T)=C(Ty—) 


Como já foi dito a determinação destas funções 
é tratada noutro local. Também noutro local é tra- 
tado o problema da determinação de K, no caso 
duma camada de ar. 


4. DETERMINAÇÃO DOS COEFICIENTES DE TRANS- 
MISSÃO DE CALOR DO AR 


4.1 Caso duma lâmina de ar 


Existem normalmente lâminas de ar no interior 
das paredes de chaminés. Estas lâminas não podem 
ser tratadas simplesmente com base em fenómenos 
de convecção. No entanto e para manter a genera- 
lidade da equação (3.6) pode definir-se uma condu- 
tibilidade equivalente. 


Devido ao facto de na generalidade dos casos a 
espessura da lâmina ser pequena em relação aos 
seus diâmetros máximo e mínimo poderemos consi- 
derar que o problema é o da transmissão de calor 
entre duas placas planas paralelas com ar de per- 
meio e teremos a geometria da figura 3. 


4 

A 

e 
E, 

/ | 


lomina de or 


a 


figura 3 


Sendo a temperatura de referência para a con- 
vecção a média de T e T, teremos que o número 
de Nusselt é dado por Nu = Ra sendo €C e n duas 
constantes dependente de Ra. 
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Ra 


2X 10:a 2,1X105 10,2 (D/H)1/9/1/4 regime laminar 


2,1X105 a 1,1X10:210,071 (D/H)1/911/3 regime turbulento 


Como Qi E h (T, — To) é-nos possível definir 


a condutibilidade equivalente K' = E 


Por outro lado a radiação contribui com um 
fluxo (ver [2] pág. 447) e supondo o ar transpa- 
rente à radiação o que é mais ou menos exacto 


4 4 
(T— To) 
iu RAE = 7 


define-se então a condutibilidade equivalente K” que 
é tal que 


e portanto 


(T4 — T$) D 
1/e, + 1/,—1 T,—T 


1 2 


K” = o 
a condutibilidade global equivalente é então 


K =K'+K” 

Os seus valores variam conforme D e a tempe- 
ratura de referência do espaço de ar. Para maior 
comodidade transcrevemos a tabela fornecida por 
[14] para tal efeito. 


Condutibilidade térmica equivalente 
[Wim K] 


Espessura da lâmina de ar (cm) 


Temperatura 
média 
(9C) 


Tabela 4.1 


4.2. Convecção do calor dos gases para a superfície 
interior da chaminé 


O número de Reynolds de transição dum escoa- 
mento interior a um tubo tem o valor aproximado 
de 2100 e a transição estende-se até aproximada- 
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mente 10000. Se notarmos que numa chaminé as 
velocidades de escoamento dos gases são raramente 
inferiores a 1Im/s e que as suas temperaturas são 
da ordem dos 300º C, somos levados a concluir que 
bastará que o diâmetro interior da chaminé seja 
superior a 10 em para que nos encontremos em regi- 
me turbulento no interior de condutas circulares. 

A mesma razão levou-nos a desprezar os efeitos 
da convecção natural por a velocidade axial ser 
suficiente para a tornar desprezável. 

Foram recolhidas & equações de Nu em função 
de Re e Pr para regime turbulento dentro duma 
conduta cilíndrica das quais 3 são teóricas. Em [5] 
obtivemos as expressões teóricas que são 


0,0396 ie” Pr 
— Prandtl Nu => ——. (4.1) 


1+2,44 Re * (Pr—1) 


que fornece bons resultados para Pr =1 


E 
0,0396 Re * Pr 

— Hoffman Nu = > 
q see 


1+1,5 pr “pe "(Pr a) 
(4.2) 


que é uma modificação da expressão de Prandtl e 
cujas propriedades são calculadas a Ta 
-— Von Karman 

0,0396 Rar Pr 
Na a api ce SE 


1+0,955 Re “f(Pr —1)+In [1 + (Pr— 17) 
(4.8) 


As expressões recolhidas que provêem de experi- 
mentação são 


— Nusselt 


0,8 0,14 1 
Nu = 0096 Re” prº (-*.) E da 


(—) (4.4) 
quo L 


0,8 
— Mac Adams Nu = 0,023 Re Pr (4.5) com n = 


0,4 fluido a aquecer 
=| 0,3 fluido a arrefecer 


que é recomendada para gases a temperaturas 
inferiores a 33º C 


— Sieder e Tate Nu = 0,027 Re” pró ( a 


us 


(4.6) 


que se pode aplicar a gases a temperatura supe- 
rior a 33ºC e para os quais 0,7 < Pr < 16,7 


— Dittus e Boelter 


0,8 0,4 
Nu = 0,0243 Re Pr ao aquecer o fluido 
(4.7) 


0,8 0,3 
Nu = 0,0265 Re Pr ao arrefecer o fluido 


(1) Ver apêndice 1. 
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aplicável quando 10: < Re < 1,2X105 


8 
Nu = 0,023 ao prº (4.8) 


— Colburn 

Para todas estas expressões calculam-se as pro- 
priedades do fluido a T, ||, exceptuando-se a equa- 
ção (4.8) para a qual apenas C, é calculado âáquela 
temperatura sendo as restantes propriedades calcu- 
ladas a T 

As expressões (4.5) e (4.6) são propostas 
por [5], enquanto [9] apresenta (4.6) e (4.4). Final- 
mente (4.7) e (4.8) são propostas por [18]. De notar 
que (4.6) é recomendado tanto por [5] como por [9] 
citando ambas o facto de Mac Adams a recomendar. 

Para verificar qual a expressão que melhor se 
pode adoptar ao problema presente escolhemos um 
caso típico em que a temperatura de homogenei- 
zação do fluido é de 200º C e a temperatura de filme 
é de 190ºC. Fizemos a representação gráfica (1) 
para valores do número de Reynolds entre 14X 10 
e 10: e depois a comparação com valores medidos 
para o ar apresentados em [5]. Daqui resultou que 
as equações que mais se parecem aproximar dos 
valores experimentais para o ar são as equações (4.6) 
e (4.7). Em função das citações recolhidas resolve- 
mos adoptar a equação de Sieder e Tate. 


4.3. Transmissão de calor num escoamento perpen- 
dicular a um cilindro 


Existem várias expressões do número de Nusselt 
em função do número de Reynolds para a trans- 
missão de calor dum cilindro para um fluido em 
escoamento perpendicular a ele. Numas os coeficien- 
tes são constantes e noutras dependem de Re, em- 
bora todas tenham sido obtidas experimentalmente. 

As seguintes soluções são citadas em [8] 


— Mac Adams 


p V D7 0,52 
Nu=LP -032+0,43 Es | 


fim Hriim 


(4.9) 


0,52 
ou Nua 0,32 + 0,43 Re 


que é válida para 0,1 < Re < 1000 e foi posterior- 
mente modificada por Eckert e Soehngen 


— Eckert e Soehngen 


p VDT7Y% 
Nu = 0,43 + 0,48 posa (4.10) 


Étiim 


que segundo [18] apresenta um desvio de 4 % em 
relação à solução de Hilpert (ver adiante). Douglas 
e Churchill afirmaram no entanto que estas duas 
correlações anteriores estavam erradas e propuseram 
uma nova correlação. Mas HSU sugere que é válida 
para Re < 500. 
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— Douglas e Churchill 


h D a AIR Ya ai, A 
Nu= 046 | | +000128, | 


film Vrilim Vfilm 
(411) 


que se aplica a regimes com 500 < Re < 109 


— Colburn 


” h Cu |% 


HS 
Prim Criim “oo K 


(4.12) 
ju encontra-se representado graficamente em [2] 
donde é possível retirar valores para interpolação 


que deixa de ser válida para Re > 105 devido a tran- 


sição de regime. 


— Ulsamer Nu = CRem Prn (413) 
Re C 
0,1 a 50 0,91 | 0,385 0,31 
50 a 10000] 0,60 | 0,500 0,31 
para a temperatura de referência T... 
— Hilpert Nu =C Rem (só para o ar) (4.14) 
Re m 
las 0,330 
4 a 40 0,385 
40 a 4000 0,466 
4 000 a 40 000 0,618 
40 000 a 400 000 0,805 


a temperatura de referência é T,,,, e estes valo- 
res foram obtidos com ar mais ou menos homogeéneo. 
Esta expressão é citada conjuntamente por [5], [9] 
e [18]. 

Para comparação foram executadas duas repre- 
sentações gráficas de (4.10), (4.11), (4.12) e (4.14), 
uma para 100 < Re < 2,5X 105 e outra, que pode ser 
considerada uma extrapolação, para 2,5X105 < Re < 
< 1,5X 108. (Ver apêndices 2 e 3). 

Destas representações podemos concluir que até 
Re = 500 qualquer das soluções é válida e que a 
partir de Re = 6000 a solução (4.10) afasta-se fran- 
camente das outras. Também para Re > 105 a solu- 
ção (4.12) se afasta das restantes como a restrição 
ao seu uso fazia prever. 

Parece, portanto, que tanto (4.11) e (4.14) são 
fiáveis. A solução de Hilpert foi por nós adoptada 
pois que interpreta mais fielmente as mudanças de 
regime que os resultados experimentais relevaram 
existir para os valores aproximados de 4, 40, 4X 103 
e 4X 104 do número de Reynolds. 
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Outras soluções como as de L. V. King [9], 
Reither [9] e uma outra de Hilpert [18] não foram 
consideradas pois os seus domínios de aplicação são 
muito restritos. 


4.4. Correcções aplicáveis às equações 


Como a temperatura não é a mesma em todos 
os pontos num mesmo instante poderá haver lugar 
para a correcção de valores do número de Nusselt 
calculados, pois as propriedades dos gases variam 
fortemente com a temperatura. Segundo Kays e 
London [10] os valores do número de Nusselt cal- 


culados a T, ,, deverão ser corrigidos por meio de 


Nu To a ua no 
E TOS FT em vez de 
Usuik - hulk Myuik 


com os valores de temperatura expressos em tempe- 
ratura absoluta e aplicável quando 


0,33 < < 3,0 


bulk 


Ainda segundo os mesmos autores, para escoa- 
mentos turbulentos completamente desenvolvidos no 
interior de tubos cilíndricos em que os gases são 
arrefecidos o expoente n toma o valor zero pelo que 
não é necessário aplicar qualquer correcção. 

Assim a equação de Sieder e Tate adoptada 
em 3.2. toma a forma 


0,8 
Nu = 0,027 Re“* pré 


Subsiste ainda o problema da influência da 
região de entrada. Das equações representadas, ape- 
nas (4.4) apresenta o factor que realiza essa correc- 
ção. Segundo Mac Adams a correcção que se aplica 
à equação de Sieder e Tate é 


h 


—— =1+ 
h 
o 


pelo que a equação de Sieder e Tate foi adoptada 
na seguinte forma final 


0,8 
Nu = 0,027 Re Pr” EA 


O efeito da região de entrada é sensível pelo que 
deve ser considerada. 


5. DETERMINAÇÃO DAS PROPRIEDADES DO AR 
Para que o estudo dos fenómenos de convecção 


no ar seja feito de uma forma precisa é necessário 
que possamos calcular as suas propriedades com 
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um certo rigor, mais a mais quando, como era o 
caso, se pretende escrever um programa para cálculo 
em computador. Assim, teve-se o cuidado de utilizar 
formas simples e facilmente calculáveis mantendo no 
entanto a precisão. 


Segundo [15] as propriedades do ar (viscosidade 
dinâmica, condutibilidade térmica e calor específico 
a pressão constante) são muito pouco influenciadas 
pela variação de pressão e os desvios que se podem 
verificar são muito pequenos para variações da 
ordem = 0,1 atm. Isto significa que os valores cal- 
culados para 760 mmHg podem ser aplicados entre 
684 mmHg e 836 mmHg. Este intervalo é muito 
maior que o intervalo de variação da pressão atmos- 
férica nos níveis baixos da atmosfera, 


Para valores de referência tomaram-se os valo- 
res tabelados em [20] para a viscosidade e para a 
condutibilidade térmica; em [13] para o calor espe- 
cífico. 


5.1. Massa específica 


Aplicou-se directamente a lei de Charles e Gay- 


-Lussac 
pT = Clte 


e de acordo com [15] Cte=, T = 1,2929 X 
X 273,15 Kg Km-s, Este valores são apresentados 
nesse trabalho tendo como fonte o Handbook of 
Chemistry and Physics, The Chemical Rubber Co. 
(1969/70). 


5.2. Viscosidade dinâmica 


E utilizada a fórmula de Sutherland sob a forma 


T,+273,15+0 
(2. | eee 
CCC Tr27315+C 


T+273,15 —/3% 
—— — —— [Ns m=<] 
T +273,15 

utilizada por [15] em que a constante C é, de acordo 
com as recomendações de V. Vasilesco, citado em 
Pascal P., «Nouveau Traité de Chimie Minérales», 
Tome I, Masson e Cie. (1956), 


= X 11 + 113 


sendo T em “C. O valor de « é 


uy = 1,716 X 105 
T = OC 


Ns m-? 


Esta fórmula provou ter desvios mínimos em relação 
aos valores tabelados em [20] e é mais exacta que 
os apresentados por [5] a págs. 37 com as constan- 


tes da tabela A. 10 da mesma e por [20] a págs 308. 
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5.3. Condutibilidade térmica 


De acordo com a teoria cinética dos gases a 
condutibilidade térmica do ar segundo [3] é 


K=e Cv qu 


Como o cálculo desta expressão não é expedito, tal 
como o processo indicado em [17], tentou-se aplicar 
a expressão fornecida em [11] que face aos valores 
tabelados em [20] apresentou desvios significativos. 
Assim resolveu-se aplicar aproximações dos valores 
de [20] por meio de polinómios. Foram ensaia- 
das as aproximações por meio de polinómios sim- 
ples pelo programa para este efeito de [8] e por 
séries de polinómios de Chebyshev segundo [4]. 
Ambas as aproximações se revelaram excelentes, 
mas devido à necessidade de se proceder a uma 
transformação de coordenadas na série de polinó- 
mios de Chebyshev preferiu-se o polinómio simples 
devido à sua maior rapidez e facilidade de cálculo 
pois o grau dos polinómios é pequeno. Assim 
K-a T'ta, T2+ta T+a, em W/mk com T 
em Kelvin e 


a — 3,2090942 X 103 
a 1,25808 X 104 
a, =— 1,1141 X 10-17 

a — 6,85634 X 10-11 


Esta aproximação é válida para o intervalo 
275 K < T < 550 K 


5.4. Calor especifico 


Na determinação do calor específico foram ini- 
cialmente testadas as fórmulas apresentadas em [7] 
e [15] cujos desvios em relação aos valores tabe- 
lados em [13] eram apreciáveis. 

Tal como para a condutibilidade térmica optou- 
-se pela aproximação polinominal que forneceu a 
seguinte fórmula 


C-aTitaT+a, Tº+a, T+a, em J/kgK 


p 
sendo 


991,610492 
0,175710352 
—  1,0434 X 103 
2,51071 X 108 
—  1,6474 X 10-9 


pop 
| 


Esta aproximação é válida para o intervalo 
275 K <T< 550K 


6. ALGUNS RESULTADOS 


Nesta secção analisam-se alguns resultados obti- 
dos pelo programa que se comparam com valores 
obtidos por outros meios. 

Foi usada como modelo uma chaminé de betão 
de acordo com [14]. 
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Caudal de gases 0,5 Kg/s 
Altura 19 m 
Diâmetro médio 0,65 m 


Espessura betão 0,15 m (base), 0,1 m (cimo) 


caixa de ar 0,05 m ( » ),0,05m( » ) 
refractário 0,05 m ( » ),005m( » ) 
Temperatura ambiente 20º C 
Temperatura dos gases à entrada 
da chaminé 200º € 
Condutibilidades térmicas betão 1W/mkK 


caixa de ar 0,6 W/mK 
refractário 0,08 W/mkK 
3 m/s 
760 mmHg 


Velocidade do vento 
Pressão atmosférica 


De acordo com os ábacos de [16] e [19] os 
gases devem sair da chaminé a 173,9ºC e 171ºC. 
Atendendo a que se trata de ábacos tomaremos um 
valor médio para referência ou seja 172,5º C. Reali- 
zado o cálculo obteve-se para a temperatura de 
saída dos gases da chaminé o valor de 184º C. Por- 
tanto a precisão que podemos obter ao utilizar 
ábacos é escassa e tem origem no facto de para o 
uso destes ábacos serem necessários apenas a tempe- 
ratura de entrada dos gases na chaminé e o diâmetro 
desta. Adiante ver-se-á que outros factores condi- 
cionam a temperatura de saída dos gases. 


6.1. Erro devido à altura das secções 


No cálculo atrás referido cada secção teve uma 
altura de 0,5m. Realizaram-se dois testes nos quais 
as secções tinham 0,2m de altura num e im noutro. 
As duas temperaturas calculadas foram iguais 
(183,9º C) e as velocidades de saída dos gases coin- 
cidiam nos três casos. 


6.2. Efeito devido à temperatura exterior 


Foi efectuado um ensaio em que a temperatura 
exterior era de 10ºC em vez de 20º C. A tempera- 
tura calculada foi 183,0ºC e a velocidade de saída 
dos gases passou de 3,886 m/s a 3,900 m/s. A tiragem 
passou de — 7,717 mmH,O para — 8,516 mmH, O. 
É possível concluir que a temperatura exterior pouco 
influencia o processo. 


6.3. Velocidade do vento 


Além do ensaio normal (a 3m/s) realizaram-se 
outros a 10m/s e a 20m/s. Os resultados podem 
comparar-se na tabela 6.1. 


velocidade do vento [m/s] 3 10 20 
tiragem [mm H, 0] —TIIT —TIO —T,ITOT 
temperatura de saída dos 


gases [ºC] 184,0 183,5 183,3 
velocidade de saída dos 
gases [m/s] 3,894 3,900 3,889 


Tabela 6.1 
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A conclusão é que a velocidade do vento pode 
ser desprezada sem que se cometa um erro apre- 
ciável. 


6.4. Influência do caudal 


Caudal [Kg/s] 0,25 0,5 0,75 1,00 3,00 
Tiragem [mm H, 0] —8,086 —7,717 —6,780 —5,395 +20,044 
Temperatura de saída 


dos gases [ºC] 169,4 184,0 189,0 191,7 197,2 
Velocidade de saída 
dos gases [m/s] 1,885 3,894 5,906 7,920 24,040 
Velocidade/caudal 7,54 7,79 7,87 7,92 8,01 
Tabela 6.2 


Verifica-se assim que uma variação de caudal pode 
trazer uma alteração bastante sensível na tempera- 
tura de saída dos gases, como aliás era de esperar. 

Os áhacos, no entanto, não estão preparados 
para este tipo de condições. 


6.5. Influência dos materiais e sua espessura 


Nos testes realizados mantiveram-se sempre o 
diâmetro interno, a altura da chaminé e as espes- 
suras do betão. Fizeram-se variar as espessuras do 
material refractário e da caixa de ar e num caso 
ainda se fez variar a condutibilidade térmica do 
refractário. 


Condutibilidade do 
refractário [W/mkK] 0,08 0,08 0,08 0,08 0,15 
Espessura do refrac- 


tário [m] 0,05 0,1 0,05 0,1 0,05 
Espessura da caixa 
de ar (m] 0,05 0,05 0,1 0,1 0,05 


Tiragem (mm H, 0] — 7,717 — 7,793 — 7,727 —7,197 —1,613 


Temperatura de saída 


dos gases [ºC] 184,0 189,8 184,8 190,1 176,3 
Velocidade de saída 
dos gases [m/s] 3,894 3,944 3,901 3,947 3,829 
Tabela 6.3 


Os resultados apresentados na tabela 6.3 provam 
que a temperatura de saída dos gases da chaminé 
está fortemente dependente da espessura do material 
refractário e da sua condutibilidade térmica, podendo 
apresentar valores bastante afastados uns dos outros. 


7. CORROSÃO EM CHAMINÉS METÁLICAS 


Na secção anterior verificou-se que a variação 
da espessura dos materiais usados (principalmente 
os refractários) e do caudal dos gases afectam sen- 
sivelmente a temperatura de saída dos gases. E 
evidente que a temperatura mínima da parede inte- 
rior da caminhé acompanha a variação da tempera- 
tura dos gases e portanto o problema torna-se crí- 
tico se se estiver a trabalhar uma gama de tempe- 
raturas próxima da zona de condensação de gases 
corrosivos. 
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Gills e Lee realizaram um estudo [21] sobre 
uma chaminé metálica com um isolamento de 1 1/27 
de lã de rocha. Este estudo compreendeu a medição 
da temperatura de saída dos gases e da temperatura 
mínima da parede interna da chaminé. Efectuaram 
também o cálculo dstas duas temperaturas por via 
teórica. 

A chaminé foi ensaiada nas seguintes condições: 


Altura útil 
Diâmetros internos 


81 ft. 
4 ft 8,5” à entrada 


[1] 


[2] 
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4 ft a 18 ft de altura 
2 ft 6º de 23 ft de altura até ao cimo 


Espessura 1.º placa de aço 2 mm 
lã de rocha 1 1/2” 
2.8 placa de aço 2 mm 
Condutibilidades aço 43,27 W/m K 
lã de rocha 0,04864 W/m K 
Velocidade média do vento 3 m/s 
Pressão atmosférica 760 mm Hg 


A temperatura dos gases à entrada da chaminé 
dependia da carga pedida à fornalha a que a chami- 
né estava ligada. Assim, os problemas de corrosão 
surgiam quando a fornalha trabalhava a carga 
reduzida. 

Os ensaios que realizámos com o programa de 
cáleulo permitiram-nos obter duas curvas de tempe- 
ratura mínima da parede interna, uma para 100 
(90ºF) e outra para -1,1ºC (30ºF). Os resultados 
obtidos levam a uma muito boa concordância entre 
os valores medidos por Gills e Lees e os valores por 
nós calculados, o que vem confirmar a grande fiabi- 
lidade dos nossos cálculos e permitem a sua aplica- 
ção prática com resultados seguros. A figura 4 é 
bem a prova disso. 

Finalmente e a título documental inclui-se uma 
tabela comparativa entre os resultados de [21] e os 
por nós obtidos. 


Plena carga 
nossos valores Gills e Lee 


temperatura à entrada 187,8ºC 187,8ºC 
Gases 

temperatura à saída 183,59C 183,9º%0 
Temperatura mínima do metal 176,29C 170,0º€ 
Carga reduzida 

temperatura à entrada 171,1%€ 1T1,IMO 
Gases 

temperatura à saída 152,0º€ 158,9+€ 
Temperatura mínima do metal 135,3%0 132,8º€ 
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Factores do dimensionamento duma chaminé 
emissora de poluentes 


RESUMO 


Neste trabalho descreve-se sucintamente o cál- 
culo da altura duma chaminé emissora de poluentes 
e acentuam-se os pontos críticos: concentrações má- 
vimas de poluentes a admitir no solo, condições 
meteorológicas e sobre-elevação da pluma. 

Como exemplo da sensibilidade do modelo às 
condições de cálculo apresenta-se o dimensionamento 
da chaminé principal da Refinaria da Petrosul em 
Sines. 


1. GENERALIDADES 


O dimensionamento duma chaminé de modo a 
assegurar a diluição adequada dos constituintes no- 
civos é um processo em que é necessário ter em 
conta um conjunto de variáveis mutuamente depen- 
dentes. 

Os principais constituintes duma emissão gasosa 
são N,. CO,, O, e vapor de água. Para além destes, 
geralmente considerados inócuos, e conforme o pro- 
cesso que dá origem à emissão da pluma, podem 
existir em pequena escala outros constituintes sus- 
ceptíveis de causar problemas a homens, animais, 
plantas ou materiais se a sua concentração exceder 
um certo valor. Um dos mais vulgares, senão o mais 
comum, é o SO, presente nos gases de combustão 
de combustíveis fósseis líquidos e sólidos. 

O SO, é o efluente cuja presença tem sido 
objecto de maior atenção não só pelo perigo que ele 
mesmo representa mas também porque a sua con- 
centração na atmosfera é normalmente tomada como 
índice do estado geral da poluição no local. 

Para que a concentração dum constituinte nor- 
malmente perigoso emitido por uma chaminé seja 
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ABSTRACT 


We briefly describe the calculation of the height 
of a pollutant emission stack. The critical points 
are stressed: maximum ground pollutant concentra- 
tion, meteorological conditions and plume rise. 

As an example of the model's sensivity to the 
input variables, we present the calculation of the 
Petrosul Refinery's main stack. 


suficientemente baixa no receptor vários processos 
se podem usar: 


1. Limitar o teor do poluente na emissão 

2. Tratar os efluentes à saída da chaminé 

3. Utilizar uma altura adequada da emissão de 
modo a aproveitar a capacidade dispersante 
da atmosfera: 
a) elevando a chaminé 


b) aumentando a velocidade e temperatura de 
saída dos gases. 


O processo mais comum e do qual nos ocupare- 
mos neste trabalho é o recurso a uma elevação 
adequada da pluma, embora em princípio a utilização 
dum processo não exclua a dos outros. 

As características de turbulência da atmosfera 
promovem a dispersão dos gases nela introduzidos. 
Quanto maior for a altura a que se dá a emissão 
de poluente, mais diluída chegará a pluma ao solo. 
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Fig. 1 


Dispersão de poluentes emitidos por uma chaminé 


Em princípio, elevando a chaminé seria sempre 
possível obter a diluíção necessária no solo. No 
entanto não será realizável uma chaminé de 600m, 
por exemplo, além de que a capacidade de dispersão 
da atmosfera diminui com a altitude e com as 
inversões de temperatura, frequentes a algumas cen- 
tenas de metros. 

Por outro lado, as muito altas chaminés (300 m), 
embora permitam evitar concentrações indesejáveis 
no solo na zona circundante, podem lançar os poluen- 
tes na circulação geral da atmosfera que os trans- 
porta a regiões distantes. Será então necessário 
limitar a emissão dos poluentes (processos 1. e 2.). 

Um caso bem conhecido é o da Suécia. Segundo 
Fisher (1975) a contribuição britânica para a depo- 
sição média anual se SO, na Suécia é 10 %, a sueca 
20 % e os restantes 70 % são provenientes de fontes 
da Europa Central e da deposição natural. 


2. MÉTODO DE CÁLCULO 
DA ALTURA DUMA CHAMINÉ 


O método de cálculo da altura duma chaminé 
envolve 3 fases, 
2.1. Estabelecimento da concentração máxima de 


poluente a admitir no solo 


Para a previsão dos efeitos nocivos duma emis- 
são gasosa tem sido comprovada a prática corrente 
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que consiste em calcular a localização e amplitude 
da concentração máxima a admitir no solo — CMS — 
de poluentes perigosos. Só em casos muito especiais 
é que esta concentração terá que ser considerada a 
nível diferente do do solo. 


A concentração observada depende do tempo 
durante o qual se faz a média, diminuindo com o 
aumento desta. 


Embora seja fundamental a existência de regu- 
lamentação sobre valores máximos de concentração 
a admitir no solo, em Portugal não existem normas 
nem sequer valores indicativos recomendados pelo 
que se tem recorrido a legislação estrangeira, 


2.2. Cálculo da elevação total da pluma 


O cálculo da elevação total da pluma — chamada 
altura efectiva da chaminé — de modo a que a con- 
centração máxima de poluente no solo tenha um 
determinado valor é efectuado através dum modelo 
de dispersão turbulenta de poluentes, considerando 


a chaminé como uma fonte pontual contínua elevada. 


Os modelos analíticos utilizados decorrem da 
resolução da equação da dispersão turbulenta intro- 
duzindo algumas drásticas simplificações mediante 
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C (x,0,0) 


as quais a equação se torna linear, obtendo-se uma 
solução analítica gaussiana (1). Os modelos mais 
correntes deste tipo são o de Sutton (1947) e o de 
Pasquill-Gifford (1961). 

Tanto num modelo como noutro a concentração 
à distância x da fonte segundo a direcção do vento 
depende de: 


1. altura efectiva da chaminé 

2. velocidade do vento 

3. parâmetros empíricos definindo a dispersão 

atmosférica nos planos horizontal e vertical 

4. taxa de emissão do poluente 

A forma geral dum modelo gaussiano que cal- 
cula a concentração no solo segundo o eixo da 
pluma emitida por uma fonte pontual contínua ele- 
vada é: 


o 
- 


qu Pla? 


va z 


C(x,0,0) = (1) 
To 


Q — taxa de emissão da pluma 

u — velocidade média do vento (m/s) 

H — altura efectiva da chaminé (m) 

9,» 0, — desvio padrão da concentração do poluen- 

temos planos transversal e vertical; são 
função da distância à fonte e estabili- 
dade aerodinâmica da atmosfera (m) 


cMs 


XKmax x 
Fig. 2 
Modelo Gaussiano — Concentração de poluente no solo, 
segundo o eixo da pluma 


Sutton apresentou a equação (1) numa forma 
diferente usando em vez dos desvios padrão da con- 
centração o, eo, os parâmetros de difusão C, eCc,e 
o parâmetro de turbulência n relacionados por 


v2 0, = C, Xx1-n/2 
VZe, = C,xt-n/2 
Os valores experimentais para C, e Cen são 


menos detalhados e completos que os de Pasquill- 


-Gifford para o, e q, que se baseiam numa série 


grande de experiências. Estes autores apresentaram 
uma classificação da estabilidade aerodinâmica da 
atmosfera em 6 classes—da extremamente instável A 
à moderadamente estável F--baseada em medidas de 
insolação e vento à superfície, sendo o, e o, apresen- 
tados sob a forma gráfica. O modelo de Pasquill- 
-Gifford tem sido extensamente utilizado em cálculos 
de dispersão de poluentes atmosféricos e de material 
radioactivo sendo recomendado pelos manuais de 
cálculo da ASME (1968), da U.S. Atomic Energy 
Commission (1968) e de Turner (1970). 

A base experimental do modelo de Pasquill- 
-Gifford é muito mais sólida e extensa do que a do 
modelo de Sutton no que se refere aos valores dos 
parâmetros de dispersão. Martin & Barber realiza- 
ram uma quantidade de trabalho experimental que 
permite relacionar os valores calculados pelo modelo 
de Sutton com os valores medidos no solo da con- 
centração máxima de SO, na região em que ocorre 
o máximo de concentração, para diferentes tempos 
de amostragem. 

Em relação ao modelo de Pasquill-Gifford não 
é do nosso conhecimento a existência de dados 
experimentais que permitam a conversão dos valores 
calculados pelo modelo para diferentes tempos de 
amostragem. Quando aplicado a situações de polui- 
ção urbana é normal considerar a concentração cal- 
culada pelo modelo de Pasquill-Gifford como a média 
de 1 hora. A directiva VDI 2289 citada por Détrie 
(1969) dá a correspondência entre as concentrações 
detectadas com tempos de amostragem diferentes 
(entre 3 minutos e 4 dias) para o SO, emitido por 
uma fonte isolada em terreno plano e sem constru- 
ções, admitindo para a concentração uma lei normal. 

Para calcular H é necessário definir as condi- 
ções meteorológicas da zona e a sua duração. Quando 
as emissões são baixas, as inversões a pequena 
altitude originam as maiores concentrações no solo 
que podem persistir mais de 1 hora. Para emissões 
elevadas, enquanto as condições de instabilidade 
atmosférica ocasionam no solo altas concentrações 
de curta duração, a atmosfera neutra com ventos 
de velocidade e direcção variáveis está na base da 
maior parte das ocorrências de níveis elevados de 
poluentes no solo. 

O cálculo da altura efectiva da chaminé H 
baseia-se na resolução da equação (1) para a con- 
centração máxima a admitir no solo e condições 
meteorológicas que conduzam às mais graves situa- 
ções de poluição. 


2.3. Cálculo da Altura Geométrica da Chaminé 


A altura atingida pelo eixo da pluma gasosa 
emitida por uma chaminé é geralmente denominada 


(1) A solução das equações completas do escoamento e dispersão turbulentos dos poluentes na atmosfera tem sido entre 
nós objecto de investigação. Foram desenvolvidos e testados métodos numéricos que permitem resolver as equações 
parabólicas às derivadas parciais que regem o fenómeno. Entre outros, citam-se os trabalhos de Domingos e Filipe 


(1972) e Santos e Filipe (1972). 
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altura efectiva da chaminé H e pode ser considerada 
como a soma de duas parcelas: a altura geométrica 
da chaminé h, e a sobreelevação 4h. 


mio consta dt Omedi sh me 


e o o 
ne o 


Fig. 3 
Sobreelevação da pluma e altura efectiva da chaminé 


H=h +ah (3) 


A sobreelevação 4h é devida a dois fenómenos 
independentes: quantidade de movimento vertical 
ligada à velocidade de saída dos gases e impulsão 
de origem térmica devido à menor densidade dos 
gases que são emitidos a temperatura superior à 
atmosférica. Conforme os casos, domina um ou 
outro fenómeno sendo comum — caso das centrais 
térmicas, por exemplo — 4h ser fundamentalmente 
determinado pela impulsão térmica. 

Os factores que afectam Ah são a velocidade 
média do vento, a temperatura da pluma, as carac- 
terísticas geométricas da chaminé e a velocidade de 
saída dos efluentes. 

A sobreelevação da pluma é gradual verifican- 
do-se, para efeitos práticos, um nivelamento a 
10 alturas de chaminé (10 h,). Como este estudo se 
baseia na estimativa da CMS, define-se altura efec- 
tiva da emissão como a altura da pluma sobre o 
ponto de concentração máxima. 

Para o cálculo de Ah têm sido desenvolvidas vá- 
rias fórmulas—cerca de 20 —-com condições de aplica- 
ção diferentes e dando muitas vezes resultados con- 
traditórios. Briggs (1969) apresenta um bom resumo 
da situação. Domingos (1973) analisa e compara 
algumas das fórmulas mais usadas que fundamentam 
os regulamentos de construção de chaminés de diver- 
sos países. No mesmo trabalho constata-se também 
que «.. examinando os vários métodos, não deixa de 
obter-se a convicção de que na maioria dos casos 
existiu subjacente uma altura de chaminé a não 
exceder e cujo valor se procurou justificar». [Na 
Fig. 4 apresenta-se a comparação das fórmulas de 
cálculo de alturas de chaminés de países diversos, 
observando-se uma dispersão apreciável dos valores 
obtidos). 

Dentre as fórmulas existentes têm conhecido 
mais sucesso as de Holland — regulamento dos 
U.S. A. (EPA) — Bosanquet, Stimke — regulamento 
alemão — Concawe, Lucas e Briggs. A fórmula de 
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Briggs — equação (4) — é baseada em análise dimen- 
sional, contrariamente à de Lucas, do Concawe e de 
Holland que são puramente empíricas. A de Briggs 
é a mais bem fundamentada teoricamente e a que 
melhor concordância apresenta com medidas experi- 
mentais para diferentes tipos de emissão. A de 
Stiimke — equação (5) — tem uma extensa validação 
experimental e apresenta bons resultados para potên- 
cias médias mas dá valores baixos para potências 
elevadas. 


ah =20Fi1i3s u-l x? (4) 
em que F= 3.8 E-5 Qu 
ah =1.5wdu-l + 65ds5/2u- (AT/T)) (5) 


àh — sobreelevação da pluma(m) 
F — fluxo de impulsão 
Qu — potência térmica emitida pela chaminé 
(cal/s) 
u-— velocidade média do vento à altura da 
emissão (m/s) 
x — distância à chaminé segundo a direcção do 
vento (m) 
d — diâmetro da chaminé (m) 
w — velocidade de saída dos gases da chaminé 
(m/s) 
AT — diferença de temperatura entre os gases à 
saída da chaminé e o ar ambiente (ºK) 
T, — temperatura dos gases à saída da cha- 
miné (ºK) 


3. APLICAÇÃO A UM CASO CONCRETO — CHAMINE 
PRINCIPAL DA REFINARIA DA PETROSUL 


O dimensionamento da chaminé principal da 
Refinaria da Petrosul em Sines é um caso com algum 
interesse e que exemplifica como se podem obter 
resultados muito diferentes a partir do mesmo mo- 
delo variando um pouco as condições de cálculo. 

Foi-nos proposto pelo Gabinete da Área de Sines 
o dimensionamento da referida chaminé e foram-nos 
facultados os dados de projecto — características 
geométricas da chaminé e da emissão e alguns dados 
sobre a meteorologia da região. 

Numa primeira fase foi acordado o cálculo da 
altura da chaminé usando o modelo de Pasquill- 
-Gifford, a fórmula de Stiimke para a sobreelevação 
e considerando a atmosfera extremamente instável — 
classe A de Pasquill. A concentração média máxima 
a admitir no solo foi estabelecida pelo GAS em 
250 ug/m: de ar para um tempo de amostragem 
correspondente ao do modelo usado sem qualquer 
correcção, i.e. 1 hora. Como seria de esperar o 
resultado foi extremamente exigente: chaminé muito 
elevada e redução drástica do teor de enxofre a 
admitir no combustível. 

Numa segunda fase repetimos o estudo em 
condições de menor instabilidade atmosférica — 
classe B--e sendo a CMS corrigida pela tabela 
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Fig. 4 
Comparação dos regulamentos de alturas de chaminé 


do VDI para um tempo de amostragem de 1 dia. 
A altura de chaminé é mais baixa e o teor de S a 
admitir no fuel-oil bastante mais elevado. 

Por outro lado foi consultada pelo GAS uma 
empresa estrangeira que, utilizando os mesmos dados 
de projecto, o modelo de Sutton, a fórmula de 
Briggs e atmosfera neutra, obteve resultados seme- 
lhantes ao 2.º estudo. 

As condições e resultados dos 3 estudos encon- 
tram-se condensados no Quadro II. 


4. CONCLUSÕES E CRÍTICA 


Os pontos críticos do dimensionamento duma 
chaminé emissora de poluentes são os seguintes: 


1. Estabelecimento das concentrações máximas 
de poluente a admitir no solo — CMS — jun- 
tamente com o tempo de amostragem e pro- 
babilidade de ocorrência. 

2. Escolha das condições meteorológicas críticas 

3. Cálculo da sobreelevação da pluma 

Como é exemplificado nos cálculos efectuados 

para o dimensionamento da chaminé da Petrosul, 
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uma alteração aparentemente pequena nos dados 
pode conduzir a variações muito importantes nos 
resultados. 

Observando o Quadro I verifica-se que com o 
mesmo modelo (Pasquill-Gifford) se obtêm resul- 
tados bastante diferentes variando um pouco a dura- 
ção do tempo de amostragem da CMS e a estabili- 
dade da atmosfera. Deste modo, mantendo-se a 
CMS de SO, pode-se quadruplicar a emissão de 
poluente desde que se considere que a estabilidade 
aerodinâmica crítica da atmosfera é representada 
pela classe B de Pasquill (atmosfera ligeiramente 
instável) em vez da classe A (estremamente instá- 
vel) e que o tempo de amostragem considerado 
aumenta a de 1 hora para 1 dia. A alteração 
geométrica da chaminé no 2.º caso é ligeiramente 
inferior à do 1.º caso (296m em vez de 300m). 

É necessário não esquecer que um erro relativa- 
mente pequeno no valor de Ah tem repercuções 
importantes no preço da construção da chaminé, 
sobretudo se esta for muito alta. No caso de uma 
altura efectiva de 4)0m se a sobreelevação for 
200m, a altura geométrica será 200m. Se o valor 
calculado de Ah for afectado de um erro de — 10 %; 
este erro traduz-se num aumento de 20m na altura 
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QUADRO | 


DADOS RESULTADOS 
SS E” o 
Classe estab. CMS Emissão 
MORRO ah atmosfera (uBg/mº) SO, (t/dia) 
1 | Pasquill-Gifford Stlimie | extremamente 250/hora 
instável 
o | Pasquill-Gifford stimke | moderadamente | oro/gja 


Sutton com coef. de Catchpole 


Briggs 


instável 


neutra 


20/dia (1) 


(1) No estudo 3. foi fixada não a concentração média mas a concentração máxima normal correspondendo a 98% de 
probabilidade de ocorrência. 


da chaminé, o que, neste caso causará um aumento 
apreciável do custo da construção. 

Finalmente queríamos acentuar o carácter apro- 
ximado dos modelos gaussianos de dispersão turbu- 
lenta dos poluentes na atmosfera e a forte compo- 
nente empírica presente nos vários passos do cálculo 
da altura da chaminé. Com efeito, a determinação 
dos parâmetros de difusão (o, o, ou C,, C,n) é expe- 
rimental assim como a determinação das constantes 
da tabela do VDI que convertem a concentração 
observada para tempos de amostragem diferentes — 
e as constantes de Catchpole que relacionam esta- 
tisticamente as concentrações calculadas pelo mo- 
delo de Sutton com as observadas para diferentes 
tempos de amostragem e probabilidades de ocorrên- 
cia diferentes. 

Cabe agora salientar que «...poderá sempre 
objectar-se que há factores de incerteza no cálculo 
que reduzem a validade das conclusões obtidas. 
A observação é certamente legítima e permite acen- 
tuar que um método de cálculo para alturas de 
chaminés deve considerar-se como de valor indica- 
tivo. A verificação «in loco» de que as concentrações 
serão excedidas deve obrigar a actuar a entidade 
oficial competente para as fazer desaparecer. A de- 
monstração de que os limites admissíveis nunca serão 
alcançados deverá permitir de igual modo a redução 
das imposições legais». In Domingos (1973). 
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Simulação numérica da congelação de peixe 
por circulação forçada de ar 


RESUMO 


Resolve-se um problema bidimensional de conge- 
lação de peixe por meio dum programa de cálculo 
numérico que permite obter a solução dum sistema 
de equações às derivadas parciais do tipo: 


com condições fronteiras do tipo 


ôu, Em fa: Rai 
= + (U ) + 0= 


5n 


pelo método ADI. 

Nesta resolução toma-se em conta que a congelação 
não se verifica a uma dada temperatura, mas sim 
numa determinada gama de temperaturas e que é 
desigualmente distribuída, conduzindo a um problema 
com propriedades físicas variando discretamente, O 
programa utilizado pode ser obtido no CTAMFUL, 
Av. Rovisco Pais, Lisboa-t. Apresentam-se também 
gráficos com a distribuição das isotérmicas que 
evidenciam o avanço da frente de congelação ao longo 
do tempo. 


t. INTRODUÇÃO 


O presente trabalho é fruto duma solicitação da 
indústria à Universidade para resolução dum proble- 
ma concreto dentro do espírito de cooperação e en- 
tendimento que se vai consolidando. Neste caso foi 
contactada a Linha de Frio Industrial do CTAMFUL 
com vista ao estudo da congelação de peixe em blocos 
num túnel de congelação por circulação forçada de 
ar. Este trabalho constitui a fase inicial de um es- 
tudo mais completo o qual inclui a realização de tes- 


(') Assistente eventual do I. S. T., Lisboa, 


J. M. BRISSON LOPES (') 


ABSTRACT 


This paper describes the numerical solution of 
a fish freezing process by means of a computer pro- 
gram which solves sets of coupled equations of the 
kind: 


du; Su, ou, 
ua ter tê 
with bondary conditions obeying to: 
du, 
ma dv (UU). L3=0 
3n + | F8 


applying an ADI method. It was assumed that free- 
2ing occurs in an interval of temperature, the rate 
of freezing varying throughout this interval of tem- 
peratures, as it would happen when freezing sea 
water. The result is a problem with discrete variation 
of physical properties which is solved for an exemple 
of fish freezing by forced air convection. The compu- 
ter program is available from CTAMFUL, Av. Ro- 
visco Pais, Lisboa-tl. 

The results are also presented in graphics show- 
ing how the freezing front and the isothermal lines 
progress throughout the freezing process. 


tes num túnel de congelação industrial. A finalidade 
do estudo é a obtenção de resultados que proporcio- 
nem aos projectistas bases para melhor compreensão 
e cálculo das instalações de congelação deste tipo. 

Com vista à realização dos testes antes mencio- 
nados foi realizada em computador uma simulação 
numérica do processo de congelação de peixe em blo- 
cos utilizando um programa de uso geral desenvolvido 
no CTAMFTUL, Este relatório pode pois também ser 


Investigador do Instituto Nacional de Investigação Científica (C. T. A. M. F. U. Lo. 
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considerado como exemplo duma aplicação a proble- 
mas concretos de programas de cálculo numérico an- 
teriormente desenvolvidos pela Linha de Métodos de 
Cálculo Numérico em Transferência de Energia e Ma- 
téria do CTAMFUL. A existência de tais programas 
gerais torna relativamente expedita a sua aplicação 
a casos particulares, sendo este o motivo por que o 
estudo preliminar do problema se iniciou pela simu- 
lação numérica do processo. 

O objectivo do trabalho é obter dados concretos 
sobre a duração do tempo de congelação, velocidade 
de congelação e propagação das isotérmicas com vis- 
ta à estimativa das necessidades posteriores de ins- 
trumentação dos testes em túnel, 


2 GENERALIDADES SOBRE O PROBLEMA 


Para a congelação em blocos o peixe é disposto 
em tabuleiros rectangulares de pequena altura que 
são colocados em estantes de vários andares sendo 
estas posteriormente colocadas no túnel de congela- 
ção. 

A simulação foi efectuada para um tabuleiro de 
tipo análogo aos utilizados nos túneis de congelação. 
Tais tabuleiros são construídos em alumínio com as 
dimensões interiores de 480 x 252x 60 mm. Tratando-se 
de paredes de alumínio, a resistência térmica da pare- 
de dos tabuleiros é desprezável face às restantes. Por 
outro lado, sendo o peixe disposto duma forma com- 
pacta admitiu-se que o mesmo forma um bloco homo- 
géneo, Com as hipóteses anteriores o problema reduz- 
-se a um problema de condução de calor em regime 
variável e propriedades físicas variáveis. O problema 
é tridimensional. Todavia a sua bidimensionalização é 
admissível se tomarmos em conta as dimensões e 
orientação dum tabuleiro face à corrente de ar frio, 
tal como a fig, 1 apresenta, 


Fig. 1 


Como o arrefecimento até ao centro do bloco é 
muito mais rápido segundo as direcções da altura e 
largura do bloco pode-se desprezar a contribuição do 
arrefecimento segundo a direcção do comprimento, 
com o que o problema se reduz a duas dimensões. 

Consideremos agora um elemento de volume in- 
finitamente pequeno. Em corpos que sofrem mudan- 
ças de fase é corrente supor que existem fontes de 
calor nos volumes elementares em que foi atingida 
a temperatura de mudança de fase considerando-se 
que a quantidade de energia libertada sob a forma de 


calor é o calor latente de mudança de fase corres- 
pondente ao volume elementar considerado. Este pres- 
suposto é válido para os casos em que a mudança de 
fase se processa a uma temperatura fixa como é o 
caso de corpos simples, No caso da congelação de 
peixe, tal já não é admissível pois trata-se da con- 
gelação de uma solução (água +sais). A solidificação 
da solução dá-se então numa gama de temperatura 
a qual está normalmente compreendida entre aproxi- 
madamente — 1ºC (este valor depende da concentra- 
cão de cloreto de sódio) e — 21,2ºC (ponto eutético 
da solução de cloreto de sódio em água). Do facto 
resultam duas consequências importantes. 

A primeira é que abaixo da temperatura do pon- 
to eutético é impossível congelar mais água e um 
abaixamento subsquente da temperatura da solução 
de cloreto de sódio dar-se-á a concentração constante 
cujo valor depende da concentração inicial da solução. 
Assim, a congelação considera-se completada quando 
a temperatura média do peixe é de —21,2ºC. 
Como do processo de congelação (por considerações 
de transmissão de calor em regime variável) resulta 
que a zona exterior do peixe se encontre sempre bas- 
tante mais fria que a interior (espinha) surgiu o 
critério prático que considera a congelação terminada 
quando se atinge — 18*€ na espinha (centro do pei- 
xe) e que a uniformização de temperatura posterior 
conduzirá a um valor próximo de — 21,2º€. 

A segunda consequência é importante para o cál- 
culo numérico: como a água contida num dado ele- 
mento de volume não congela totalmente à mesma 
temperatura, não se poderá supor que as proprieda- 
des físicas se mantém constantes até ao momento da 
congelação e durante o arrefecimento subsequente a 
esta. 

Também não é possível distribuir o calor latente 
de solidificação de maneira uniforme por uma maior 
ou menor gama de temperaturas. A sua distribuição 
dependerá sim da taxa de formação de gelo que é 
função da concentração inicial da solução e do anda- 
mento da linha líquida no diagrama de fases da solu- 
ção. Assim as propriedades físicas neste intervalo de 
temperatura são função da quantidade de gelo for- 
mado e da respectiva taxa de formação. Este proble- 
ma será retomado na secção 4. 

Outra hipótese deste trabalho diz respeito à difu- 
são da solução através dos tecidos devido a diferen- 
ças de pressão osmótica, Observações realizadas por 
Plank (ver [6]) mostram que o fenómeno é negligen- 
ciável desde que a velocidade de avanço da frente de 
congelação seja suficientemente rápida. Esta hipótese 
será adoptada neste trabalho preliminar. 


3 PREPARAÇÃO DA SOLUÇÃO NUMÉRICA (ADIMEN- 
SIONALIZAÇÃO) 


31 — Equação da difusão do calor 


O programa utilizado [1] resolve numericamente 
sistemas de equações parabólicas não lineares em três 
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variáveis independentes, sendo uma delas o tempo. 
As equações assim resolúveis são forma geral 


du, õu, Su, 
EO, — Te uriçõ- JOS ——me— pm s (1) 
5t 1 5xi i2 5vi i3 
em que 
Su 
ij =— (ota, ça ) 
o e 
com ael(x,y)e i=l,..,mn 
j=1, 2,83 
K=1, cc D 


admitindo condições fronteiras do tipo 


Yu) O (2) 
E 


sendo 5, y e é funções de x, y e t. 


As equações devem ser previamente adimensio- 
nalizadas. Para tal escolheram-se as variáveis adi- 
mensionalizadas seguintes 


Y 
— espaço SE esp oii e 8 = ET 


— tempo 


A equação da difusão do calor em coordenadas 
cartesianas a duas dimensões espaciais e na ausência 
de fontes de calor toma a forma 


5T | ST “a 37 
EC — = k—— )- ES 


Quando as propriedades físicas são variáveis (ver 
[3]) a forma da equação que nos interessa é 


JT To! 
REL ( | Rana 


St 5xº DA a 


5k sr ) ST Nº 4 

JLC PMES ii 

É conhecido que a adimensionalização nos con- 
duz aos seguintes resultados 


ST LT XL 4 BT 
— ESDACO E gem So a, . ——- , — 
E Ldxe Co Bo Pas 
Í i Í i 
3T 3T 
— tempo A ua aco 
St E dz 


pelo que (4) Será 


5T dk 35T NY 5T 
É q pm, À asaço , 
Copa = ee + gr [aee ) + (ye ) | 
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ou sob uma forma de acordo com a expressão ge- 
ral (1) 


3T LL faTANO 
crua Yu |! 5) 
Gi ( dx”, e q. | 


O cálculo das derivadas de 1.º ordem não causa 
transtorno especial pois o programa tem subrotinas 
para tal efeito. 

O cálculo dos valores da difusibilidade e da con- 
dutibilidade térmicas será devidamente considerado 
na secção 4. 


32 — Condições aos limites 


Não nos deteremos sobre a condição inicial (tem- 
peratura no início do arrefecimento) pois é simples 
informar o programa a seu respeito. Uma subrotina 
especial lê a distribuição inicial de temperatura pelos 
pontos da malha rectangular aqui considerada. Esta 
distribuição é uniforme no nosso caso mas o progra- 
ma admite qualquer distribuição inicial de tempera- 
tura. 


Para melhor compreensão das condições frontei- 
ras detenhamo-nos na fig. 2 que representa uma faixa 
dum tabuleiro na direcção do escoamento de ar à qual 
se sobrepôs um sistema de eixos 


Fig. 2 


E evidente que se as condições fronteiras forem 
do mesmo tipo em todos os pontos da fronteira é 
possível, devido à simetria, separar cada um das qua- 
tro partes. 

Assim teremos as seguintes condições fronteiras 

— fluxo nulo para x*=0o ou y*=o 

ou seja B=y=4=o em (2) 
— fluxo imposto para x*=1 ou y*=a 


o que implica que 


5T 
-=Kk— =h[(Tijs — Ta] (6) 
an 
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que é uma condição de acordo com (2). Para que se 
mantenha a semelhança térmica entre o fenómeno 
real e o modelo adimensional há que fazer com que 
Nu=Nu* o que nos impõe que se defina um coefi- 
ciente adimensional h* tal que 


adimensionalizando também em relação à normal, 
será 


ST 
a — — h8 [(* es | 
K ns h*[( Ts [a 


ou seja 
E a a E hºTa = 
ne o no pl anái 
ou ainda A=K 
p= 
= — h* Ta 


Está pois definido o problema. Resta agora cal- 
cular as propriedades físicas, o que faremos seguida- 
mente. 


4. PROPRIEDADES FÍSICAS 


Pelo que já foi dito, as propriedades físicas do 
peixe estão condicionadas pela quantidade de gelo 
formado. É este o primeiro problema a ser tratado 
visto que dele depende a determinação de todas as 
propriedades físicas a partir do início da congelação. 


41 — Gelo formado e sua taxa de formação 


A quantidade de gelo que se encontra formado a 
uma dada temperatura compreendida entre OC e 
— 21,2*%C depende da concentração inicial da solução 
de cloreto de sódio em água (visto que a influência 
dos outros sais é desprezável na água do mar) e do 
andamento num diagrama de fases da linha líquida 
de tal solução, 

A bibliografia consultada sobre este assunto ([2], 
[4], [5], [6] e [7]) fornece-nos muitos elementos, 
umas vezes sob a forma de tabelas e gráficos, outras 
de expressões. Entre as expressões encontradas há 
uma devido a Riedel que mostrou fornecer valores 
médios dos elementos confrontados. Essa expressão 
estabelece uma relação entre a fracção de água que 
está congelada a uma dada temperatura e essa mes- 
ma temperatura, é 
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| O quando T> 272,15 K 


1,105 
ti) + | 0,31 
“og (27315 — T) 


se 2722 45K>T >251,95K 


0.900 se T< 251,95 K 


A taxa de formação de gelo é então 


0 se TH 272,15K ou TL251,9 
de 0,34255 
dT (27345 — T) In 10 [log (273,15 — T) + 0,31 


se 272 145K > T > 251,95K 


42 — Calor específico 


Como o gelo deixa de formar-se a temperatura 
constante, a equação da difusão do calor não pode 
ser resolvida com a introdução de fontes de calor a 
uma dada temperatura correspondente ao ponto de 
congelação. A equação terá que ter a forma (4) e 
assim somos levados a definir um calor específico 
equivalente em que se tenha em conta o gelo formado 
e a respectiva taxa de formação, o que nos leva a 


de 


dT 


C= Co(1-—-W)LCaW(1— 0) + CgoWo—ILNW 


que pode ainda escrever-se sob a forma 


des a j 
O - Cs = W Co —-— Ca nem E, RS 4 to Cs si Ca )|] 
dT 


onde a formação de gelo é responsável pelo termo 


au. 
sã dT 


43 — Condutibilidade térmica 


A condutibilidade térmica terá uma expressão 
que também resulta da ponderação das propriedades 


K=KU-W+KWA-D+K, Wo 
ou 
K=K,+W [K—K,+o (K, — K,)] 


donde resulta 


dK dá 
SA wik—Ki 
dT (Re— Ka) 


44 -— Valores numéricos adoptados 


Para a simulação numérica que realizámos esco- 
lhemos um modelo formado por peixe gordo tal que: 
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W = 0,6 


—3 —3 
p = 980 k, m  =0,98 g em 


Para as restantes propriedades os valores adoptados 
foram 


L =7 9,91 kcal kg! 
Cs = 0,440 kcal kg 
Ca = 1,009 kcal kg) 
Cg = 0,460 kcal kg' 
Ks = 0,22 kcal him? Kº = 
Ka = 0,487 kcal him? K+ = 
Kg = 1,905 kcal him? K* = 


= 334521><10M erg g" 
K+ = 1,8 44><107 erg g* Ka 
K+= 42244812 erg g'! K+ 
K!= 19259280 erg g" K+ 
25586 erg stem? Ki 
56 597,2 erg s'cmi Kº 
221 542,4 erg stcm* Kº 


5. RESOLUÇÃO NUMÉRICA 
5.1 — Condições e adaptação ao modelo 


Trata-se dum modelo adimensionalizado consti- 
tuído por uma malha rectangular igualmente espaça- 
da nas duas direcções do espaço. Segundo uma das 
direcções, x, o modelo é adimensionalizado à dimen- 
são 1. Como os tabuleiros cuja congelação pretende- 
mos simular têm as dimensões internas 480X 251 X 
X 60 mm, a simulação irá realizar-se para um quarto 
de secção cujas dimensões são 126X 30 mm. Se optar- 
mos por 6 nós segundo x teremos uma distância 
de 0,2 entre dois nós consecutivos. Assim ter-se-ão 
22 nós na direcção y. 

A temperatura de evaporação da instalação onde 
se vão realizar os testes é de —40ºC, e existe normal- 
mente uma diferença de temperatura de 5ºC entre o 
ar e o fluido refrigerante, Podemos então supor que 
a temperatura não perturbada do ar é de —35€. 
A temperatura do peixe a congelar à entrada do túnel 
pode ser suposta uniforme e igual a 15 €. 

Para coeficiente de transmissão de calor por con- 
vecção (que neste caso é forçada) tomámos um valor 


—1 º —2 — 
aproximado de 22,7 kcal h m K para o seu 


valor médio. Assim 


o o | —1 —2 ..—l 
h=22,7kcalh m K <=26388,75ergs cem K 


5.2 — Resultados 


Os resultados que fundamentalmente nos inte- 
ressavam eram a velocidade de congelação, os tempos 
de congelação e de refrigeração e a distribuição de 
temperaturas ao longo do processo. 

Assim e dos resultados obtidos que incluiam a 
distribuição de temperaturas com intervalos de apro- 
ximadamente 2 minutos, foi possível obter os seguin- 
tes valores 


Tempo de congelação ( —18%C no centro) 
3h 3min. 425 

Tempo de congelação ( -—21,2X de média) 
2h  Sâ2min. 25s 
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Tempo de refrigeração ( —1º%C de média) 
45 min. 10s 
Início da congelação na menor espessura 
4min. 47s 
Atinge-se —590 no centro ao fim de 
2h  3émin. 52s 
Velocidade de congelação 
1,2 em/h 


Como se vê a velocidade de congelação corres- 
ponde à de uma congelação ainda lenta, mas já pró- 
xima do limite a partir do qual é considerada rápida, 
o que nos permite concluir que a hipótese feita sobre 
a difusão de água foi correcta. 


Com base nas distribuições de temperatura que 
o programa calculou traçaram-se as curvas isotérmi- 
caz para instantes com aproximadamente meia hora 
de intervalo entre si, Esta traçagem foi realizada com 
um programa traçador de curvas de nível desenvol- 
vido no CTAMFUL para o calculador HP9830A uti- 
lizado no CTAMFUL o qual está acoplado a um 
traçador de gráficos. 

Os gráficos obtidos evidenciam bem o avanço da 
frente de congelação e das isotérmicas e permitem 
uma boa visão do seu andamento e distribuição. 
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SIMBOLOGIA T, — Temperatura não perturbada do fluido circulan- 


te (ar) 
Letras T, —'Temperatura inicial do peixe 
u, — Variável genérica da equação diferencial i 

C — Calor específico x — Variável espacial 
C, — Calor específico da água y — Variável espacial 
C, — Calor específico do gelo W — Fracção de água na massa total do peixe 
C, — Calor específico dos tecidos secos do peixe 
h — Coeficiente de transmissão de calor por conven- 

ção Letras gregas 
K — Condutibilidade térmica 
K, — Condutibilidade térmica do gelo a — Difusibilidade térmica 
K, — Condutibilidade térmica dos tecidos secos do a. — Difusibilidade térmica inicial 

peixe p — Massa específica do peixe 
L — Calor latente de fusão da água Tr —'Tempo adimensional 
| — Dimensão característica do modelo real, factor o — Fracção de água congelada a uma dada tempe- 

de escala ratura 
n — Normal à fronteira do sistema dirigido para o 

exterior 
Nu — Número de Nusselt Indice superior 
T -—'Temperatura absoluta 
t — “Tempo real * — Variável adimensionalizada 


Fig. 3 
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Fig. 4 


Isotérmicas às Oh 30miln. 47,55 
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Isotérmicas às O h 59 min. 31,9 s 
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Fig. 6 
Isotérmicas à 1h 30min. 19,45 
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Fig. 7 
Isotérmicas à 1h 59 min. 3,6 s 
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Fig. 8 
Isotérmicas às 2h 36 min. 0,2 5 


Fig. 9 
Isotérmicas às 3h 4min, 44,2 s 
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As células fotovoltaicas no aproveitamento da energia solar 


RESUMO 


A energia solar representa uma solução que se 
considera possível para resolver em parte o grave 
problema de encontrar partes de energia que possam 
substituir em termos de futuro as disponibilidades 
energéticas clássicas que num futuro mais ou menos 
próximo se encontrarão esgotadas. 

Neste artigo tecem-se algumas considerações 
acerca do aproveitamento da energia solar, usando 
células solares as quais convertem, a radiação solar 
directamente em energia eléctrica com uma razoável 
eficiência. 


1. INTRODUÇÃO 


Tomando como factor reconhecido a importância 
que obrigatoriamente se tem de dar ao aspecto de 
encontrar novas formas de energia para assegurar 
em termos de futuro as disponibilidades de fontes de 
energia, num planeta como o nosso em que as reser- 
vas clássicas são limitadas e as necessidades sempre 
crescentes, pode e deve por-se a questão de qual o 
caminho a seguir, quais as formas alternativas de 
energia a explorar dentro do que se apresenta como 
possível. 


Uma das hipóteses é sem dúvida o aproveita- 
mento da energia solar, dadas as potencialidades que 
na realidade apresenta, em países como o nosso. 


Sem nos debruçarmos muito sobre o pormenor 
das vantagens do aproveitamento da energia solar, 
podemos no entanto referir que a quantidade de ener- 
gia que atinge a superfície da Terra em cada ano é 
de cerca de 3X 10* Joules, número que convém com- 
parar com as actuais reservas orgânicas conhecidas, 
que se estimam em 1X1024 J, com um consumo pre- 
visto para o ano 2000, de cerca de 69xX1010 G Joules. 
E este o momento do nosso país se debruçar seria- 
mente sobre o problema, pois uma perda de tempo 
neste campo, traduzir-se-á num irrecuperável atraso 
relativamente a outros países que preparam já a 


LEOPOLDO J. M. GUIMARÃES (') 


SUMMARY 


Solar energy represents one possible solution to 
the major problem we face today: how to find energy 
resources to replace the non-renewable energy reser- 
ves we have been using and wasting up to now. 

In this paper some considerations are made on 
solar cells which are photovoltaic cells that convert 
sunlight directly into electrical energy with a good 
efficiency. 


industrialização da fabricação de centrais solares para 
exportarem para países considerados subdesenvolvi 
dos, os quais, por sua vez, usufruem normalmente 
das condições climatéricas mais favoráveis a este tipo 
de aproveitamento. 

A propósito, indicamos a seguir os níveis de in- 
vestimento de alguns países na investigação e desen- 
volvimento da energia solar, sob vários aspectos do 
seu aproveitamento: 


Milhões de $U.S. 


Estados Unidos TO 
Japão 36 
França o 6 
Alemanha Ocidental 6 


2. ALGUMAS CONSIDERAÇÕES SOBRE CÉLULAS 
SOLARES 


A célula solar é uma célula fotovoltaica que con- 
verte energia da radiação solar directamente em 


(!') Professor Extraordinário da Universidade Nova de Lisboa, 
Investigador do Instituto Nacional de Investigação Científica (C. T. A. M. F. U. L.). 
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energia eléctrica, com uma eficiência que à luz do 
actual desenvolvimento tecnológico atinge valores não 
muito superiores a 10%*. A potência eléctrica gerada 
é aproximadamente proporcional à intensidade da ra- 
diação solar e à área útil colectora, não exigindo a 
conversão qualquer movimento mecânico, nem ciclos 
termodinâmicos intermédios. 

Fabricando-se as células solares a partir de ma- 
teriais semicondutores, interessa referir alguns aspec- 
tos que influenciam a obtenção de resultados de inte- 
resse prático. 

Duma forma simplista, os materiais semicondu- 
tores podem ser classificados a partir da sua largura 
de banda, parâmetro que permite conhecer o rendi- 
mento energético teórico, se forem feitas algumas 
aproximações. A Fig. 1 resulta dum estudo feito nesta 


Rena. 


Teórico 


10 (é (8 2.2 
Fig t (ev) 


base, podendo ver-se que certos materiais se encon- 
tram junto do máximo da curva de variação. 


Os rendimentos obtidos na prática, divergem mui- 
tos dos valores teóricos, principalmente por razões 
tecnológicas. Cite-se o exemplo do GaAs (junção p-n 
obtida por difusão num monocristal) que apresenta 
um rendimento de apenas 11%. 


Existem de facto vários factores que limitam a 
eficiência duma célula fotovoltaica. Uns dependem 
da própria fenomenologia do funcionamento da célula, 
de modo que pouco se pode esperar no sentido de 
melhorar os seus efeitos; outros são passíveis de se 
tornarem sucessivamente menos importantes, à medi- 
da que se progride na técnica de fabricação. 

De entre os primeiros, podemos destacar: 


— Utilização duma parte muito restrita da ener- 
gia dos fotões incidentes responsáveis pela 
criação dos pares electrão-lacuna, sendo a 
energia excedente perdida através do meca- 
nismo de criação de fonões. 


— Absorção incompleta da radiação incidente. 
— Factor de tensão, dado pela relação entre a 
tensão do curto-circuito e a diferença de po- 


tencial correspondente à descontinuidade ener- 
gética característica da junção p-n. 


26 Largura de banda 


De entre os segundos, salientamos: 


— Perdas por reflexão óptica na superfície. A mi- 
nimização deste efeito pode ser conseguida 
através da aplicação de camadas anti-reflec- 
toras sobre a superfície colectora, Tem-se aliás 
avançado bastante neste campo, com imediatas 
repercussões na fabricação de módulos solares 
que fazem o aproveitamento da energia solar 
transformando-a em energia térmica. 


— Recombinação de pares electrão-lacuna ou seja 
incompleto aproveitamento de pares que difun- 
dindo-se através da junção, não contribuem 
para o mecanismo de conversão de energia. 


— Radiação perdida devido à pequena espessura 
da parte útil da célula, com consequente absor- 
ção da radiação através dos contactos exte- 
riores. 

— Degradação das características da junção, de- 
vido à resistência interna. 


3. MATERIAIS PARA FABRICAÇÃO DE CELULAS 
SOLARES 


As primeiras células utilizadas, manufacturaram- 
-se a partir de monocristais de Si e pouco depois a 
partir de cristais de CdS. Muito embora o Si permita 
células estáveis com eficiências superiores a 10%, 
grandes facilidades na obtenção da matéria-prima a 
par duma tecnologia bastante mais evoluída que a de 
qualquer outro material, a sua vulnerabilidade à dete- 
rioração pela radiação e a exigência de monocristais 
de elevada pureza, portanto muito caros, fazem com 
que se dispense uma atenção sempre crescente a ou- 
tros materiais que apresentam potencialmente melho- 
res condições de utilização. Concretamente não se 
conseguiu até agora substituir com vantagem global 
a tecnologia do Silício, mas está naturalmente aceite 
que os novos materiais apontam como hipóteses pro- 
missoras para células solares de segunda geração. 
Estão neste caso materiais como CdTe, Cus, InP, 
AsGa que correspondem normalmente a materiais se- 
micondutores caracterizados por terem uma largura 
de banda que cai entre 1 ev e 1,6 ev. Esta largura de 
banda, corresponde uma optimização da potência for- 
necida pela célula, ou seja a um valor máximo do 
produto tensão disponível pela corrente gerada. 

As células solares podem apresentar-se em forma 
monocristalina ou em forma de películas finas. Muito 
embora a eficiência destas últimas, seja bastante infe- 
rior às primeiras em face do menor tempo de vida 
dos transportadores minoritários, apresentam no en- 
tanto algumas vantagens que as consagram como 
possíveis soluções para o futuro. Algumas dessas ca- 
racterísticas são as seguintes: 


(*) A eficiência teórica máxima que se pode conseguir para uma célula solar fixa-se em 25% em condições de radiação 
directa da luz solar num material de largura de banda de 1,4 ev. 
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— maior área da célula 

— menor peso 

— maior resistência à degradação por efeito de 
radiação 

— menor preço dos conversores, quer para uso 
terrestre quer espacial. 


4. CONSIDERAÇÕES FINAIS 


As células solares tiveram a sua principal apli- 
cação no fornecimento de energia a módulos espaciais, 
dadas as suas características de reduzido peso, sim- 
plicidade, eficiência e facilidade de utilização em 
qualquer parte do espaço. Note-se, porém, que a 
curto prazo, se prevê a sua utilização em larga escala 
em aplicações terrestres, podendo mesmo afirmar-se 
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que a breve trecho, digamos já na década 1980-1990, 
a associação de todas as formas de aproveitamento 
de energia solar, poderia cobrir pelo menos 20% das 
nossas necessidades em energia o que evidentemente 
torna indispensável conjugar esforços para estabele- 
cer ao nosso nível um programa coerente de inves- 
tigação. Existem de facto pessoas e grupos a tentar 
trabalhar neste domínio mas duma forma perfeita- 
mente desgarrada sem uma planificação adequada 
que permita retirar o melhor rendimento do esforço 
investido. Naturalmente que esta situação resulta em 
nossa opinião, da forma lacunar como se tem enca- 
rado a investigação aplicada no nosso país, não se 
estabelecendo as prioridades necessárias por um lado 
e não se dotando essa mesma investigação das ne- 
cessárias verbas que constituiriam sem dúvida um 
investimento que obrigatoriamente produziria a seu 
tempo os necessários frutos, 
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Elaboração de tabelas de temperatura Sol-Ar para paredes, 
para a região de Portugal Continental 


SUMÁRIO 


De todos os métodos utilizados para o cálculo 
de cargas térmicas em Ar Condicionado, o da Tem- 
peratura Sol-ar é o mais preciso e é largamente 
utilizado. Por estes motivos e dado que as tabelas 
existentes não se referem à situação climatológica 
de Portugal, elaborou o autor, tabelas de tempera- 
tura Sol-ar próprias para as características do 
nosso país. 


INTRODUÇÃO 


Em qualquer projecto de ar condicionado, um 
dos mais importantes problemas que se põem na 
sua resolução é o cálculo das cargas térmicas a que 
está sujeito o ambiente a climatizar, ou seja a 
energia que esse espaço está a trocar com o exterior. 
A sua importância é imediata já que o dimensiona- 
mento de toda a aparelhagem é feita com base no 
valor calculado da carga térmica. Ter-se-ã portanto 
que, se a carga térmica calculada for inferior à que 
realmente existiria para as condições pretendidas, 
a instalação terá no final uma potência inferior 
âquela que é realmente necessária. Por outro lado, 
se a carga calculada for muito superior à carga 
real máxima (para as condições pretendidas), em- 
bora a instalação possua a potência necessária, a 
sua zona de funcionamento, muito abaixo da do 
rendimento óptimo, torna muito difícil ou até impos- 
sível um perfeito controlo das condições interiores 
dentro dos limites pretendidos. 

Se, por um lado, um método simples tem a 
vantagem dum cálculo rápido da carga térmica, a 
precisão do resultado é muito difícil de prever à 
priori, já que não se entra em conta com as proprie- 
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SUMMARY 


The Sol-air temperature method is an accurate 
and widely used method in Air-Conditoning engineer- 
ing, for the calculation of heat gains in a room, 
This work presents tabulated values that are to be 
used in Portugal, as evaluated tables of Sol-air 
temperature for this weather conditions are not 
available. 


dades dos materiais que constituem a estrutura do 
espaço a climatizar. Um grande desvio em algumas 
destas propriedades em relação às que serviram de 
base às fórmulas empíricas que constituem o algo- 
ritmo de cálculo, podem levar a uma das situações 
acima mencionadas, com o inconveniente de o reco- 
nhecimecato de tal erro só muitas vezes poder ser 
detectado depois de montada toda a instalação, 


Pelas razões indicadas, dos métodos conhecidos 
pelo autor, o da temperatura Sol-ar é o mais preciso 
e portanto o aconselhável, visto ser aquele que entra 
em linha de conta com o maior número de proprie- 
dades dos materiais de construção através dos quais 
se farão sentir os fenómenos de transmissão de 
calor. No entanto não é do seu conhecimento que 
existam tabelas próprias para o clima de Portugal 
e é no sentido de tabelar os valores da temperatura 
Sol-ar para as condições do nosso país, que foi 
desenvolvido este trabalho, tomando como base para 
os valores da radiação solar, os apresentados por 
P. Moon em Proposed Standard Solar Radiation 
Curves for Engineering Use. 
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MÉTODO DA TEMPERATURA SOL-AR 


Para melhor compreender este método, supo- 
nhamos a parede P do edifício representado na 
figura 1. As trocas de calor entre o ambiente (A) 
e o exterior (E), devem-se essencialmente a um 
gradiante de temperatura entre o exterior e o inte- 
rior e à radiação solar que incide na parede. 


Esta é composta por: 


R, — Radiação solar directa 

R, — Radiação solar difusa 

R, — Radiação solar reflectida (por objectos 
próximos) 


Supondo o edifício isolado, não existirão por- 
tanto corpos próximos e a radiação solar reflectida 
será proveniente apenas da radiação reflectida pelo 
solo. 

Embora a radiação solar extra-terrestre (radia- 
ção solar antes de atingir a atmosfera terrestre) 
seja praticamente constante, a intensidade dessa 
radiação ao atingir a superfície da Terra variará 
com a altura do Sol, já que o caminho percorrido 
pelos raios solares será tanto menor quanto mais 
próximo ele estiver do zénite (ver a figura 2). 
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Desta forma, a absorção de parte da energia 
solar (pelos gases raros, vapor de água, ozono, etc.) 
será tanto maior quanto menor for a altura solar. 
Na figura 3 está representada esta variação da inten- 
sidade da radiação solar com a altura solar. 


Intensidade da 
radiação (N/m?) 3 


00 


DD 2 “4 


70 Bo 90 
alfara solar (º) 


to 20 30 40 SO 60 


Fig. 3 


1 — Radiação solar directa e difusa. 2 — Radiação solar 

directa. 3 e 4 — Radiação solar difusa. Curvas 1, 2 e 3 de 

acordo com «Proposed Standar Solar Radiation Curves for 
Engineering Uses, 


Curva 4 considerada como sendo igual a 11.5% da radia- 
cão solar directa. 


A radiação difusa resulta de parte da energia 
radiante ser difundida em todas as direcções, ao 
encontrar moléculas de azoto, oxigénio, vapor de 
água, poeiras, gotas de água (núvens), etc. Tam- 
bém a sua intensidade varia com a altura solar e 
a proporção entre ela e a radiação solar total não 
será constante. No entanto a variação é pequena, 
como se pode vir na figura 4. 


Ra 
x, 
Mp + dm 
ff 
to 
va/or medio 
9 
to 2Zo 30 40 so 60 70 Bo go 
altura solar(?! 


Fig. 4 


Em média o seu valor corresponde a aproxima- 
damente 10,0% da radiação solar total (considerando 
apenas a radiação directa e difusa). Isto corresponde 
a poder-se considerar que a radiação difusa tem um 
valor menor ou igual, incidindo numa parede vertical 
a 5,75% da radiação solar directa ao longo do dia. 
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Aliás a intensidade da radiação difusa só é superior 
a valor apresentado e desde que se admita um erro 
de 5 % no cálculo do seu valor, para alturas solares 
inferiores a 7º ou superiores a 76º, De notar que em 
Portugal a altura solar nunca excede 76.3º, já que 
a latitude mais baixa do território é cerca de 37.2º 
e para as zonas do norte, do centro e grande parte do 
sul do país a altura solar ao longo de todo o ano não 
excede os referidos 76º. Em relação aos valores de al- 
tura solar inferiores a 7º é conveniente referir que 
essa situação apenas se verifica durante cerca de duas 
horas ao longo de todo o dia: a primeira até uma 
hora depois do nascer do Sol, e a segunda uma 
hora antes do pôr do Sol, e como já foi referido, a 
estas horas a radiação solar total é bastante pequena. 


De salientar também que não existe qualquer 
problema em considerar para cálculos, um valor 
da energia recebida pela parede por radiação difusa 
ligeiramente superior ao valor real pois não se esta- 
ria em nenhum dos dois inconveniente apontados 
no início deste trabalho. 


Quanto à radiação incidente no solo que em 
parte será reflectida, também ela depende da altura 
solar, já que a intensidade da radiação que incide 
numa superfície horizontal é tanto maior quanto 
maior for a altura solar. A radiação reflectida pelo 
solo que é a parte que nos interessa considerar 
varia também com o clima da região pois a factor 
de reflexão dum solo depende, para além do próprio 


tipo de solo, das condições climatéricas existentes 
no local. 


Convém ainda referir que uma parede por estar 
a uma dada temperatura, emite, ela também, radia- 
ção, mas esta é praticamente compensada pela radia- 
ção que essa parede recebe, emitida pelo solo vizinho. 


Pelo que foi dito é evidente que um cálculo da 
carga térmica devida à radiação solar, a partir dos 
valores da intensidade da radiação solar total e das 
diferentes orientações que os raios solares vão tendo 
em relação à parede em estudo, obrigava a cálculos 
relativamente moroso que não se coadunam com um 
método de cálculo de cargas térmicas que se pre- 
tende tão expedito quanto possível. O método da 
temperatura Sol-ar cria então um aumento fictício 
na temperatura do ar exterior e a carga térmica 
devida à radiação é considerada como uma carga 
também ela devida a uma diferença de temperatura. 
A carga térmica através duma parede, devido à 
radiação e à diferença de temperaturas entre o 
exterior e o interior será então igual à carga tér- 
mica que existiria não havendo radiação solar mas 
estando o ar exterior a uma temperatura superior 
à real. Esta temperatura fictícia do ar é a tempe- 
ratura Solar (T.). 


A diferença entre a temperatura Sol-ar e a 


temperatura exterior, que é portanto proporcional à 
carga térmica recebida pela parede por radiação, foi 


São estes diferentes valores e para diferentes meses 
que serão apresentados neste trabalho. 


CÁLCULO DA TEMPERATURA SOL-AR 


A intensidade da radiação solar directa que 
incide numa parede é função da altura solar (w) e 


q — latitude; 4 — declinação; Az — azimute; a — altura so- 
lar; H — ângulo horário; PN — Pólo Norte; PS — Pólo Sul; 
Z — Zénite; N-—- Nadir; N—- Norte; S— Sul; E — Este; 
W — Oeste; HH'— plano do horizonte; EE'— plano do 
Equador; Sar declinação solar máxima, correspondente 
ao solstício de Verão. Para este dia, o Sol nasce no ponto 1, 
ao meio-dia situa-se na posição 2 e põe-se no ponto 3. 
Para uma posição qualquer do Sol (S) estão índicadas as 
principais coordenadas tanto equatoriais como horizontais, 
e que definem a posição do Sol. 
A latitude do território do continente situa-se entre 37.2º 
e 42º. A declinação solar varia entre 23.5º N (Verão) e 
23.5º S (Inverno). 


do ângulo azimutal da parede (Azp). O ângulo 
azimutal da parede é o ângulo formado pelo azimute 
solar e pela orientação da parede. Por analogia 
com o sentido seguido para o azimute (ver figura 5) 
a orientação das paredes far-se-á também no sentido 
dos ponteiros do relógio e com origem no Norte. 
Assim, a orientação duma parede norte será 0º, duma 
parede este 90º, duma parede sul 180º, etc. Utilizando 
esta forma de proceder, o valor do ângulo azimutal 
duma parede é obtido directamente por uma sub- 
tracção de dois ângulos: azimute solar e orientação 


calculada para diversas situações ao longo do dia. da parede. 
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A partir da figura 6 obtém-se: 


I = I cosacos Azp 


sendo I, a intensidade da radiação que incide na pa- 
rede é I a intensidade da radiação solar. 


Fig. 6 
Uma vez conhecida a radiação recebida pela 
parede é possível determinar o acréscimo fictício de 


temperatura T.. 


Ata = E ai 17h, (a lg —elIç) 


sendo: 


e — coeficiente de emissividade da superfície 
exterior da parede 


290) 


a — coeficiente de absorção da superfície exte- 
rior da parede 


h, — coeficiente de convecção exterior 


I. — radiação recebida pela superfície exterior 
da parede 
IL, — radiação emitida pela superfície exterior da 


parede 


Como já foi referido, pode-se considerar 1, = 0. 


Para a elaboração das tabelas considerou-se uma 
velocidade do vento igual a 12Km/h e que o clima 
da região é um clima temperado. Estes valores são 
normais para o Verão na quase totalidade do terri- 
tório continental. 


A cor da superfície exterior da parede, para os 
valores apresentados, é a cor escura, O acréscimo vir- 
tual de temperaturas assim calculado, para diferen- 
tes alturas solares e ângulos azimutais da parede é 
apresentado na tabela 1. 


Notar que as horas referidas são horas solares. 
Aliás qualquer referência a uma dada hora será 
sempre feita em relação à hora solar e não à hora 
legal. 


Se as paredes em estudo tiverem uma cor clara 
ou média, deve-se aplicar uma correcção aos valores 
apresentados, multiplicando-os respectivamente por 
0.65 ou 0.85. Para os edifícios de cor muito clara 
usar também o factor de correcção 0.65, já que ao 
fim de algum tempo o escurecimento da parede 
provocará uma maior absorção da energia radiante 
que sobre ela incide, situando-se o seu coeficiente 
de absorção, neste caso, próximo do valor próprio de 
uma parede de cor clara. 
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TABELA 1 


SS 


ALTURA 


SOLAR 


30.6 130.5 


30.3 |30.3 


ÂNGULO AZIMUTAL DA PAREDE 


A — Parede em sombra e solo próximo ao Sol. 


AA — Parede em sombra e solo próximo em sombra. 


Se a velocidade do vento (v) for diferente da 
indicada, deve-se proceder também a uma correcção: 


. 64—v 
AT Sd ST a 


sendo a velocidade do vento expressa em Km/h e 
àT'.. O valor lido na tabela. 

Se a altitude do local for superior a 300m, a 
correcção a aplicar é: 


aT = AT” + 0.00023 (a — 300) AT” 
sa sa sa 


sendo a altitude do lugar (a) expressa em metros. 

De notar que para o caso de paredes em sombra 
estão assinaladas duas situações: a primeira refe- 
re-se ao caso em que o solo próximo está todo ao 
Sol; no segundo caso o solo está em sombra. Numa 
situação real se não houver objectos próximos poder- 


TÉCNICA 437 


13.8 [14.1 


—emm 
— 0. 


14.9 114.2 


13.5 /12.9 


12.0 [11.6 


10.5 110.3 


6.1 | 65 É 
1.9 | 2.0 | 2.0 


-se-á dizer que o acréscimo fictício de temperatura 
(AT,.) duma parede à sombra está situado entre os 
dois valores indicados na tabela, já que parte do 
terreno vizinho estará ao Sol e outra parte à sombra. 
Como é lógico a diferença de valores apresentados 
na tabela deve-se apenas à incidência na parede da 
radiação reflectida pelo solo. A vantagem desta 
forma de proceder em relação à apresentação de 
um único valor para paredes em sombra, é mani- 
festa, sobretudo nos instantes que se seguem à 
parede ficar em sombra ou antes de começar a 
receber radiação directa. Em qualquer destas oca- 
siões o valor AT da parede é superior ao que seria 
apresentado numa tabela geral. 

No entanto esta tabela, não torna ainda expedito 
o cálculo da carga térmica para uma dada parede. 
Para o seu cálculo, a uma dada hora, torna-se neces- 
sário conhecer a temperatura Sol-ar média da parede 
o que implica a necessidade de fazer um cálculo do 
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valor da temperatura Sol-ar para todas as horas, 
desde o nascer ao pôr do Sol. 

Para tornar mais rápida a utilização deste mé- 
todo, foram calculados para os meses de Março a 
Setembro os diferentes valores de AT e para as 
principais orientações de paredes. Apenas são indi- 
cados estes meses pois no território nacional, o mês 
de projecto de Verão corresponde sempre a um deles. 
Geralmente o mês de projecto em Portugal é Julho 
ou Agosto. 

Para a elaboração das tabelas, considerou-se 
uma latitude de 40º N para o cálculo da altura e azi- 
mute solares ao longo do dia, e portanto, para a 
determinação da intensidade da radiação solar inci- 
dente. Embora a declinação varie com o dia do ano, 
ela foi considerada constante ao longo de cada mês, 
variando apenas de mês para mês. O seu valor 
corresponde à declinação característica do mês que 
se está a considerar, e que é, para os meses tabe- 
lados, a existente nos dias 22 de Março, 20 de Abril, 
21 de Maio, 21 de Junho, 23 de Julho, 24 de Agosto 
e 22 de Setembro. 


Foi já tomado em conta o efeito atrás referido 
duma parede em sobra que está a receber a radiação 
solar reflectida dum solo que estará normalmente 
uma parte ao Sol e outra à sombra. 

A última linha referente a cada parede, indica 
o valor médio diário de AT. para cada uma das 
orientações apresentadas. 

A coluna referente à parede em sombra diz 
espeito a uma parede que nunca receba radiação 
directa, bem como o solo que lhe está próximo. 
Notar que esta situação é diferente da de uma 
parede orientada a Norte e que usualmente se costu- 
ma designar por «à sombra». Neste caso ela pode 
receber radiação solar directa: a seguir ao nascer 
do Sol e até cerca das 6 horas, e antes do pôr do 
Sol desde as 18 horas aproximadamente. Em rela- 
ção ao terreno que lhe está próximo, haverá sobre 
ele incidência solar em maior ou menor percentagem 
durante todo o dia, excepção feita para as 12 horas, 
hora a que estará totalmente em sombra. 

As correcções a aplicar a estas tabelas são aná- 
logas às indicadas para a tabela 1, no que se refere 
à velocidade do vento, à cor e à altitude. 


TABELA 2-A 
MES DE JUNHO 


Parede N Ss SW W NW Sombra 
0500 | 41 | 78 | 73| 29 | 06| 05 | 05| 06 | 04. 
0600 8.6 “ Bl ass IgA go | do 
0700 8.2 3.7 24 | 18 | 31 | 19. 
"0800 | 57 5.2 3.5 1.9 | 36 | 20. 
0900 | 61 “103 | 49 | 24 | 36 | 21 
"1000 | 62. “22 | 17] 63| 33| 31| 21. 
a 1100 | 47 s 175 | 86 | 53 | 28 | 21. 
é | 1200 2.1 185 54 81 1 Eu 
m | 1300 47 17.5 | 208 | 166 q | 21 
1400 | 62. “147 | 22 | 237 | 138 | 210 
1500 | 61 | 36 | 24 | 49 “103 | 249 | 287 | 195 | 21. 
1600 | 517 5 | 52 | 230 | 307 | 235 | 20. 
1700 8.2 “37 | 184 | 287 | 24 | 19. 
1800 — "86 | 20. “21 [118 | 213 | 202 | 16 
1900 | 41. 05 | 06 | 29 | 73 | 78| 04. 
Média | 37. “58 | 7o | 76 | 58 | 11 
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TABELA 2-B 
MES DE MAIO OU JULHO 


Parede N NE wW NW Sombra 
"19 | 39 | 38 | 16 | 03| o2| 03| o03| 02. 
“180 | 195. 4 | 15 | 3 
"66 | 231 | 282 | 189 | 35 | 24 | 17 | 28 | 16. 
“223. 1.9 33 | 1.8 
“183 | 288 | 261 | 117 | 50 | 25 | 33 | 20 
"8 | 125 | 239 | 254 | 163 | 64 | 35 | 28 | 21. 
ma. 5.6 25 | 21. 
5.0 8.0 50 | 21 
2.5 "221 | 166 | 74 | 21 
“1400 | 58 | 28 | 35 | 64 | 163 | 254 | 239 | 125 | 21. 
“s7 | 33| 25| 50 | 17 | 261| 288 | 183 | 20. 
3.3 | 306 | 223 | 18 
2.8 "282 | 231 | 16 . 
“mo | 18 | 12) 13| 19 | 109 | 195 | 180 | 11. 
“19 | 03| o03| o2| o03| 16| 38| 39 | 02. 
“sa | 58) 73] 731 5 “q3 | 53 | 10 


TABELA 2-C 


MÊS DE ABRIL OU DE AGOSTO 


Parede 
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TABELA 2-D 
MÊS DE MARÇO OU DE SETEMBRO 


Parede 
0700 
0800 
— mo” 
1000 | 
7 1100 25 | 151. 
o 1200 28.8 | 222 
" 1300 | 275 | 276 
1400 | 36 | “21 | 303 | 
1500 | 18.6 | 299. 
1600 | us 


Quanto à correcção em latitude, ela não é neces- 
sária. Para qualquer latitude do território conti- 
nental, os valores de AT | não diferem em mais 
de 5% dos valores apresentados para a latitude 
de 40º, excepção feita para alguns casos situados 
nas horas próximas do nascer ou do pôr do Sol, ou 
seja quando a diferença entre a temperaturas Sol-ar 
e a temperatura exterior é muito pequena, Além disso 
a diferença entre os valores de AT | para a latitude 
de 40º e para a mesma hora, não tem sempre o 
mesmo sinal, para as diferentes orientações das pare- 
des. Nestes casos ter-se-á portanto que para uma 
sala com quatro paredes exteriores, o erro cometido 
no cálculo da carga térmica, pelo método da tem- 
peratura Sol-ar é inferior a 5%, já que os erros 
cometidos pela leitura de AT de cada uma das 
paredes tenderão a anular-se. (Isto, como é evidente, 
supondo que não são cometidos quaiquer outros erros 
no cálculo da carga térmica). 

Notar que o erro de 5% foi o erro admitido 
para o cálculo da diferença de temperaturas Sol-ar, 
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11.9 | 25,5 


o que torna irrelevante a diferença entre os valores 
apresentados na tabela 1, e os que existiriam para 
outra latitude qualquer do território continental. 

A título informativo, são apresentados nas tabe- 
las 3A e Bos valores de ST, para os meses de Março, 
Junho e Setembro, meses dos equinócios e do solstício 
de Verão, para as latitudes máxima e mínima do 
território nacional, isto é, 42º N e 37.2º N. 

De referir ainda que todos os programas neces- 
sários à elaboração das tabelas foram escritos em 
FORTRAN IV tendo sido utilizado, para a sua reso- 
lução, o computador IBM 360 no sistema de explo- 
ração DOS. 
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TABELA 3-A 
MES DE JUNHO 


Lat. 37.2º Lat. 42.0º 
Parede N E s W N E S W 
0500 | 29 | 52| 04| o04| | 49 | 87| o07| 06. 
"0600 | 85 | 204 | 20 | 13 “87 | 219 2.2 1.4 
otroo | 87 | 284 | 36 | 18 | | 78 | 289 38 | 18. 
"0800 | 66 | 306 52 | 19. “b2 | 807 | 5% 1.9 
0900 | 62 | 288 | 91 | 23) | 59 | 286 | w1 | 25. 
1000 | 65 | 239 | 133 | 30 | | 59 | 236 | 157 | 34 
« | 1100 53 | 167 | 160 | 48 | | 44 | 165 | 186 | 55. 
o | 1200 | 21 | 82 | 169 | 82) “21 | 80 | 197 | 81 
* [1300 | 53 | 48 | 160 | 167 | 44 | 55 | 186 | 165 
"1400 | 65 | 30 | 133 | 239 “| 59 | 34 | 157 | 236 
1500 | 6.2 23 | 91 | 288 | - 59 | 25 | 111 | 286 
"1600 | 66 | 19 | 52| 306) | 52 | 19 | 57304 
1700 | 87 | 18 | 36 | 284) | 78 | 18| 38 | 289 
"1800 | 85 1.3 20 | 204 8.7 14 | 22 | 29 
1900 | 29 | 04 | 04 | 52 49 | 06 | 07| 87 
Média 3.8 7.4 4.8 al í|Is7rl mal se| gg 
TABELA 3-B 
MÊS DE MARÇO OU SETEMBRO 
Lat. 37.2º Lat. 42.0º 
Parede N E Ss W E 8 W 
14 | mó! 47l ããil ÕõÃITagZg |O! ATI SO. 
32 | 268 | 117 1.9 “30 | 258 | 120 | 18. 
"269 | 182 | 27. "36 | 260 | 189 | 27. 
"226 | 234 | 40 “| 35 | 29 | 245 | 40. 
155 | 267 | 56 | | 28 | 149 | 280 | 53. 
“esl 29 | col |zol c4l|zos| 65. 
“56 | 207 | 155 | | 28 | 53| 280 | 149 
“40 | 234 | 226 “| 35 | 40 | 245 | 29. 
27 | 182 | 269 | - "36 | 27 | 189 | 260. 
19 | 117 | 268 “30 | 18 | 120 | 258. 
1,1 7 | 190 | | 17 | 10 | 47 | 180. 
14 | 55 | 82 | 55. “13 | 58| 86 5.3 
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